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ABSTRACT 


Electronic  support  measures  (ESM)  systems  play  an  increasingly  important  role  in 
modem  warfare  and  can  influence  the  outcome  of  a  military  engagement.  The  application 
of  ESM  can  be  extended  to  anti>guernlla  and  anti-dmg  operations  where  law  enforcement 
agencies  can  exploit  the  fact  that  their  presence  is  inducing  the  outlaw  to  depend  more  on 
radio  communications  to  coordinate  their  activities.  When  a  propagation  path  of  no  more 
than  one  reflection  at  the  ionosphere  (l*hop)  can  be  assumed,  position  of  an  HF  emitter  can 
be  determined  by  a  single  observing  site  using  vertical  triangulation,  provided  that  the  height 
of  the  ionosphere  at  the  point  where  the  radio  wave  is  reflected,  can  be  determined.  This 
teclmique  is  known  as  "high  frequency  direction  finding  single-site-location"  (HFDF  SSL). 
This  thesis  analyzes  the  HFDF  SSL  error  in  measuring  the  direction  of  arrival  ot  the  signal, 
how  this  error  is  generated  by  the  antenna  array  and  its  effect  on  emitter  location.  The 
characteristics  of  the  two  antenna  arrays  used  by  a  specific  HFDF  SSL  system  that 
implements  the  phase-interferometer  technique  were  studied  using  electromagnetic  modeling. 
Results  showed  that  angle-of-arrival  errors  for  the  high  band  array  were  less  than  0.5®  and 
were  under  0.2®  for  the  low  band  antenna  system  The  maximum  HFDF  SSL  lateral  and 
range  error  of  this  system  were  found  to  be  8.7  Kin  and  22.4  km  respectively  for  the  high 
band  array,  when  the  targeted  emitter  is  located  at  500  kni,  and  the  incoming  wave  ha.s  an 
E-propagation  mode.  The  smallest  lateral  and  range  error  were  found  to  be  2.1  km  and  4.3 
km  respectively  for  the  low  band  array,  when  the  targeted  emitter  is  located  at  300  km,  and 
the  incoming  wave  has  an  E-mode. 
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I.  INTRODUCTION 


A.  OVERVIEW 

Passive  electronic  warfare  (EW)  systems  play  an  increasingly  important  role  in 
modem  warfare.  Electronic  support  measures  (ESMi  are  now  used  to  accomplish  a 
number  of  missions  that  include  the  detection,  location,  identification  and  threat 
assessment  of  hostile  emitters.  Modem  ESM  systems  provide  critically  important  data 
which  can  influence  the  outcome  of  a  military  engagement. 

The  application  of  ESM  can  be  extended  to  anti-guerrilla  and  anti-drug  operations. 
Law  enforcement  agencies  can  use  these  techniques  to  help  locate  illegal  installations  and 
to  obtain  additional  advantage  from  the  element  of  surprise.  Specifically,  law  enforcement 
agencies  can  exploit  the  fact  that  their  presence  is  inducing  increased  mobility  of  the 
outlaw,  which  in  turn  will  cause  them  to  depend  more  on  radio  communications  to 
coordinate  their  activities. 

Passive  ESM  systems  are  harder  to  build,  and  accuracy  requirements  are  generally 
more  stringent  than  requirements  for  active  EW  systems.  These  systems  have  to  work  in 
an  unfriendly  and  totally  uncooperative  signal  environment.  In  order  to  be  effective  they 
must  exhibit  broad  spectral  coverage,  high  probability  of  intercept,  real  time  operation, 
and  high  accuracy.  Probably  the  most  important  component  of  modem  ESM  systems  is 
the  radio  direction  finder. 
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Radio  direction  finding  (DF)  can  be  defined  as  geographically  fixing  the  position 
of  a  targeted  transmitter  both  in  azimuth  and  range,  using  one  or  more  receiving  sites. 
Virtually  all  DF  systems  derive  emitter  location  from  the  measured  arrival  angle  of  the 
received  signal.  This  angle  of  arrival  is  always  defined  in  azimuth  and  in  some  cases 
elevation.  Horizontal  or  azimuth  triangulation  is  the  determination  of  the  position  of  an 
emitter  using  azimuth  angles  measured  at  two  or  more  DF  stations.  This  is  the  most 
common  technique  used  in  DF,  but  requires  that  receiving  sites  be  separated  by  a  distance 
comparable  to  the  distance  to  the  emitter,  and  their  operation  must  be  precisely 
synchronized  in  order  to  obtain  accurate  position  measurements. 

In  the  high  frequency  (HF)  band,  propagation  of  radio  waves  takes  place  mainly  by 
means  of  reflections  at  the  ionized  layers  of  the  atmosphere.  Only  at  very  short  distances 
are  the  direct  ray  and  surface  waves  of  importance.  By  using  this  special  characteristic, 
another  way  of  DF  can  be  implemented  for  HF  signals  provided  that  the  height  of  the 
ionosphere  at  the  point  where  the  radio  wave  is  reflected,  can  be  determined. 

When  a  propagation  path  of  no  more  than  one  reflection  at  the  ionosphere  (1-hop) 
can  be  assumed  for  a  received  HF  signal,  emitter  position  location  can  be  determined  by 
a  single  observing  site  using  vertical  triangulation.  This  technique  is  known  as  "high 
frequency  direction  finding  single-site-location"  (HFDF  SSL)  and  is  used  to  determine 
location  of  HF  emitters  up  to  approximately  2,000  km  where  the  sky  wave  signal  may 
arrive  via  a  1-hop  propagation  path. 

In  the  process  of  measuring  the  direction  of  arrival  of  a  signal  using  an  HFDF  SSL, 
two  basic  errors  can  occur,  one  due  to  the  uncertainty  of  measuring  the  exact  direction 
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of  arrival  of  the  signal,  and  the  other  due  to  path  deviations  which  make  the  direction  of 
arrival  different  from  the  true  direction  of  the  transmitter.  The  first  error  depends  on  the 
accuTiicy  of  techniques  and  equipment  used,  ana  the  second  on  the  precision  in 
detemiming  characteristics  of  the  ionosphere  at  the  rcflecdon  point,  at  time  of  measure* 
ment. 

This  thesis  wUl  analyze  the  HFDF  SSL  error  in  measuring  the  direction  of  arrival 
of  the  signal,  how  this  error  is  generated  by  the  antenna  array  and  its  effect  on  emitter 
location.  The  characteristics  of  the  two  antenna  arrays  used  by  a  specific  HFDF  SSL 
system  that  implements  the  phase-interferometer  technique  will  be  studied  to  find  their 
effects  on  the  HFDF  SSL  enor. 

B.  DF  TECHNIQUES 

Direction  finding  systems  can  be  classified  in  two  broad  categories  according  to  the 
manner  in  which  their  antenna  systems  obtain  and  process  infonnation  from  the  arriving 
signal  of  interest.  These  categories  are  scalar  DF  and  vector  (or  phasor)  DF  systems 
[Pef.  1]. 

Scalar  DF  systems  measure  only  one  component  of  the  signal  of  interest  (such  as 
amplitude  or  phase),  and  use  the  information  of  the  measured  component  to  determine 
azimuth  and/or  elevation  angles  of  arrival  of  the  signal.  Vector  DF  systems  measure 
simultaneously  two  or  more  components  of  the  received  signal  and  try  to  determine  the 
characteristics  of  the  individual  rays  that  compose  the  incoming  wave-field.  A 
disadvantage  of  this  kind  of  system  is  that  a  single  incident  wave  may  be  defined  by  as 
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many  as  four  parameters  [Ref.  1];  thus  it  must  have  a  sufficient  number  of  antenna  ports 
and  sufficient  measuring  and  processing  ci^abilities  so  that  the  desired  number  of 
unknown  parameters  can  be  resolved. 

Scalar  DF  systems  are  generally  much  simpler  than  vector  systems  and  therefore 
are  better  suited  for  tactical  use.  Current  scalar  DF  systems  can  be  classified  into  three 
basic  types: 

1)  Spinning  DF 

2)  Doppler  DF 

3)  Interferometer  DF 

1.  Spinning  DF 

These  systems,  also  called  rotating  DF  systems,  measure  the  amplitude  of  a 
signal  received  by  a  highly  directive  antenna  or  antenna  array.  The  receive  pattern  of  the 
antenna  array  is  rotated  mechanically  or  electronically  until  a  null  or  peak  is  directed 
towards  the  incoming  signal.  The  angular  position  of  the  antenna  (or  its  pattern)  provides 
infoimation  to  compute  the  line  of  bearing  to  the  transmitter.  Spinning  DF  systems  are 
by  far  the  most  commonly  used  type  of  DF  systems. 

Most  spinning  DF  systems  use  a  small  aperture  antenna  (small  antenna  size 
compared  to  a  wavelength)  in  order  to  achieve  a  receive  pattern  with  sharp  nulls,  but  this 
technique  has  the  disadvantage  of  a  low  overall  gain  and  the  null  occurs  where  the  gain 
is  at  minimum.  To  overcome  this  limitation  the  patterns  of  tv/o  or  more  elements  can  be 
combined  to  produce  sharper  beamwidths  and  allow  the  use  of  the  maximum  gain  portion 
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of  the  antenna  receiving  pattern.  To  achieve  this,  small  aperture  antenna  elements  must 
be  used  and  the  overall  gain  is  still  low. 

Even  though  large  aperture  antenna  arrays  have  been  used  to  improve  the  gain 
of  the  system,  spinning  DF  systems  still  have  the  disadvantage  of  not  being  instantaneous 
since  the  amount  of  time  the  antenna  beam  points  to  a  certain  azimuth  is  a  small  fraction 
of  one  complete  rotation.  This  characteristic  of  spinning  DF  can  induce  errors  in  short- 
duration  signals  if  the  scan  rate  is  not  high  enough.  For  a  high  scan  rate  the  system  has 
to  have  a  wide  bandwidth. 

2.  Doppler  DF 

Doppler  DF  systems  constitute  a  variation  of  the  spinning  DF  technique  where 
frequency  modulation  is  superimposed  on  the  received  signal  by  mechanically  rotating  a 
single  omnidirectional  antenna  element  rapidly  around  the  perimeter  of  a  circle  [Ref.  2]. 
The  frequency  deviation  or  Doppler  Frequency  is  proportional  to  the  scan  rate  and  the 
size  of  the  circle  of  rotation. 

A  common  version  of  this  technique  obtains  an  electronically  simulated 
Doppler  effect  by  smoothly  commutating  between  fixed  identical  antenna  elements 
positioned  around  a  ring.  The  requirement  to  rotate  (or  commutate)  the  antenna  elements 
results  in  the  disadvantage  of  not  being  instantaneous. 

An  additional  requirement  is  that  the  diameter  of  the  antenna  circle  must  be 
several  wavelengths  in  order  to  obtain  adequate  resolution.  This  characteristic  makes 
Doppler  systems  better  suited  for  very  high  frequency  (VHF)  and  above. 
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3.  Interferometer  DF 


DF  systems  based  on  interfeiometeis  compare  values  of  a  specified  signal 
parameter  measured  on  two  identical  antennas  (or  groups  of  antennas)  to  derive  the  angle 
of  arrival  (AOA)  of  the  signal.  The  most  common  parameters  used  for  Interferometer  DF 
measurements  are  amplitude,  time  of  arrival,  and  phase. 

DF  Interferometers  use  the  fact  that  the  measurement  of  the  selected  parameter 
at  each  antenna  would  yield  the  same  value  if  the  signal  arrived  equally  spaced  from  both 
antennas,  creating  a  difference  null  for  this  direction  of  arrival  of  the  signal.  Any  other 
condition  will  yield  a  difference  output  value  which  is  proportional  to  its  divergence  from 
the  null  case. 

Even  though  the  three  parameters  mentioned  are  intimately  related,  the  most 
common  interferometer  used  for  DF  is  the  phase-interferometer.  This  study  will 
concentrate  on  a  specific  implementation  of  this  technique  on  an  HFDF  SSL  system.  The 
characteristics  of  a  phase-interferometer  will  be  covered  in  greater  depth  in  the  next 
chapter. 

C.  PURPOSE  OF  THE  THESIS 

The  purpose  of  this  thesis  is  to  investigate  the  magnitude  of  the  error  induced  in  a 
phase-interferometer  HFDF  SSL  system  by  the  antenna.  The  research  will  include  the 
analysis  of  two  HFDF  SSL  antenna  arrays  and  the  DF  algorithm  employed  by  the 
Southwest  Research  Institute  of  San  Antonio,  Texas  (SWRI)  in  its  land-based  HFDF  SSL 
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system.  The  Numerical  Electromagnetics  Code  (NEC)  will  be  used  to  model  the  antenna 
characteristics,  array  geometry,  and  ground  conditions  and  to  simulate  real  DF  operation. 

To  provide  framework  for  the  analysis  of  the  system,  a  tactical  scenario  will  be 
developed.  This  scenario  will  include  the  characteristics  of  medium  range  military  land 
operations,  as  well  as  land  based  anti-guerrilla  and  anti-drug  operations. 

Specific  areas  of  investigation  described  in  this  paper  include:  (1)  discussion  and 
analysis  of  phase-interferometers  and  the  effects  of  ionospheric  conditions  on  HFDF  SSL 
errors;  (2)  analysis  of  the  characteristics  of  the  antennas  selected;  (3)  modeling  and 
analysis  of  the  antenna  array  geometries  used  by  SWRI;  and  (4)  analysis  and  comparison 
of  theoretical  and  simulated  outputs  of  the  HFDF  SSL  system  and  the  determination  of 
the  magnitude  of  the  errors  generated  by  the  antenna  arrays  and  their  effect  on  emitter 
location. 

D.  TACTICAL  SCENARIO 

The  following  theoretical  scenario  was  developed  for  this  research.  When  possible, 
it  includes  the  characteristics  of  typical  medium-range  military  land  operations,  as  well 
as  land-based  anti-guerrilla  and  anti-drug  operations.  The  area  of  operations  is  located 
near  the  equatorial  region  and  the  terrain  conditions  include  the  highlands  (with  altitudes 
that  range  from  6,000  to  10,000  feet)  and  the  Amazon  jungle.  The  system  is  intended 
to  be  used  all  year  round  and  any  time  of  day  or  night. 

The  area  over  which  the  system  (including  the  antenna  array)  is  to  be  deployed  will 
usually  be  constrained  to  a  square  of  less  than  300  by  300  meters.  The  DF  antenna  array 
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should  allow  for  easy  assembly  and  disassembly,  so  it  can  be  transported  to  other 
locations  when  required. 

The  expected  distance  to  target  emitters  will  range  from  300  to  500  km  and  the 
operating  frequencies  will  range  from  3  to  15  MHz.  Most  of  the  target  emitters  will  be 
ground-based,  but  occasionally  air-based  emitters  may  be  of  interest.  The  type  of  signals 
the  DF  system  is  expected  to  work  with  are  single  sideband  (SSB)  and  frequency-shift 
keying  (FSK). 


8 


II.  HIGH  FREQUENCY  SINGLE-SITE-LOCATION  DIRECTION  WNDING 


A.  THE  CONCEPT  OF  A  PHASE-INTERFEROMETER 

A  phase-inteiferometer  is  a  DF  system  that  uses  an  airay  made  up  of  one  pair  (or 
more)  of  identical  antenna  elements  and  determines  the  direction  of  arrival  of  a  received 
signal  from  the  relative  phase  measured  between  two  antenna  elements.  Tlie  antennas  used 
in  this  kind  of  systems  are  generally  omnidirectional  or  have  broad  beamwidths  (on  the 
order  of  90  degrees  or  more).  The  difference  in  the  phase  received  by  two  antennas  which 
are  equally  displaced  from  the  incoming  wave  will  be  zero,  creating  a  phase  difference 
null  along  the  boresight  of  the  interferometer  antenna  pair. 

In  a  phase-interferometer  DF  the  signal  is  effectively  received  at  the  phase  center 
of  each  antenna  and  the  position  of  the  phase  center  relative  to  each  antenna  is  a  function 
of  the  angle  of  arrival  of  the  wave  front.  If  the  anteimas  are  identical  the  distance  between 
the  phase  centers  and  the  slope  of  the  line  connecting  them  is  independent  of  the  angle 
of  arrival  of  the  wave  front. 

The  line  connecting  the  phase  centers  of  two  antennas  in  a  phase-interferometer  DF 
is  called  a  baseline.  An  interferometer  can  have  more  than  one  baseline.  An  interferometer 
system  that  employs  only  2  antennas  is  called  a  single  baseline  interferometer  and  if  it 
employs  more  than  2  anterma  elements  it  is  called  multiple  baseline  interferometer. 

Interferometers  have  the  advantage  of  not  requiring  pattern  scanning  as  the  previous 
DF  systems.  They  are  particularly  useful  for  short-duration  signals. 
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Before  the  study  of  details  of  a  phase-interferometer  the  concept  of  a  "plane  wave" 
must  first  be  defined.  Since  electromagnetic  waves  propagate  as  an  expanding  sphere  as 
they  move  away  from  the  transmitter,  the  shape  of  the  wave  front  describes  a  circle  in 
any  plane  that  contains  the  transmitter.  At  long  distances  (many  wavelengths  away),  the 
case  of  interest  for  HFDF  SSL,  this  wave  front  can  be  treated  as  a  straight  line.  When  it 
reaches  two  antennas  whose  separation  is  many  times  smaller  than  the  distance  to  the 
emitter,  it  will  be  referred  to  as  a  plane  wave. 

Baseline  orientation  is  defined  as  the  angle  between  the  baseline  and  geographic 
true  North.  For  convenience  and  to  simplify  the  analysis,  true  North  is  considered  to  be 
aligned  with  the  positive  "x"  axis  of  a  3-dimensional  coordinate  system,  as  shown  in 
Fig.  1. 
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B.  PHASE-INTERFEROMETER  BASELINE  TYPES 

i 

1.  Single  Baseline 

A  phase-interferometer  DF  system  that  en^loys  only  two  antennas  is  called 
a  single  baseline  interferometer  and  its  geometry  is  shown  in  Fig.  2.  In  this  interferometer 
configuration,  the  incident  electric  field  intensity  at  antenna  1  will  be 

£,  =  £  ,  (i) 


where 

E  =  magnitude  of  the  electric  field  intensity, 

P  =  free  space  propagation  constant  -  In! %  and 
X  =  distance  traveled  by  the  wave. 

Similarly  the  electric  field  intensity  at  antenna  2  will  be 

£2  =  E  ,  (2) 


where 

d  =  distance  from  antenna  1  to  antenna  2  and 
\jr  ;=  elevation  angle  of  arrival  of  the  plane  wave. 

Only  the  relative  values  of  phase  are  of  interest  so  antenna  1  can  be  taken  to 
be  the  reference  and  the  electric  field  intensities  become  £,  =  £  and  £,  =  £ 
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a 


PLANK  WAVE 


Therefore 

_fl_  _  g^COrt  (3) 

and  the  signal  received  at  antenna  1  leads  the  signal  received  at  antenna  2  by  the  phase 
angle 

<I>,j  =  pdcosO  .  (4) 

From  this  relation  the  angle  of  arrival  of  the  plane  wave  can  be  calculated  by 

e  =  .  (5) 


By  making  the  length  of  the  interferometer  baseline  larger,  the  angle  0  can  be 
determined  with  greater  precision  and  the  effect  of  phase  measurement  errors  on  <^,2  wUl 


2 


decrease.  On  the  other  hand»  a  greater  number  of  ambiguities  is  introduced  since  can 


only  be  measured  within  2n  radians.  For  0  <  A  0,2  <  2n, 


and 


COS0  s 


A  0,2  +  N{2%) 


d  A  O,, 

1  (COS0)  = _ li  +  N  . 

X  27C 


(6) 


(7) 


The  second  term  of  this  last  expression  is  called  the  "interferometer 
ambiguity",  since  it  gives  the  number  of  ambiguities  that  are  introduced  in  the  single 
baseline  interferometer  by  the  choice  of  baseline  length,  d/X.  From  Fig.  3,  for  no 
ambiguities  to  occur,  a  baseline  length  of  less  than  one>half  wavelength  is  needed. 


Fig.  3.-  Phase  Interferometer  Ambiguity  Function  Diagram. 
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Thus  in  a  single  baseline  interferometer,  an  ambiguity  of  ±7t/2  about  the 
baseline  occurs  because  of  symmetry  when  the  antenna  phase  centers  are  separated  by 
more  than  one-half  wavelength.  Additional  ambiguities  occur  for  elements  spaced  at 
distances  larger  than  one  wavelength  because  the  value  of  6  starts  repeating  when  the 
phase  difference,  d>,j,  exceeds  180®. 

To  resolve  the  second  kind  of  .unbiguity  the  antenna  phase  centers  are  required 
to  be  at  a  distance  less  than  one-half  wavelength.  For  the  first  ambiguity  to  be  resolved 
the  field  of  view  of  the  antenna  array  can  be  limited  by  using  subsidiary  front  and  rear 
pointing  antennas  and  comparing  their  amplitudes.  Alternatively,  a  second  baseline  can 
be  used  and  their  phase  differences  combined. 

In  a  single  baseline  interferometer  when  the  elevation  angle  of  the  incident 
wave  is  not  zero,  the  locus  of  points  that  produce  the  same  phase  delay  forms  a  cone  of 
revolution  al  at  the  axis  of  symmetry  between  antennas  1  and  2.  The  exact  direction  of 
arrival  of  the  incoming  signal  cannot  be  determined  in  this  case  because  various 
combinations  of  azimuth,  <|> ,  and  elevation,  \{r,  can  give  the  same  value  of  phase  delay. 

Therefore,  a  two-element  phase  system  provides  the  possibility  of  measuring 
the  angle  of  arrival,  B ,  of  an  elevated  wave,  but  in  an  ambiguous  way.  There  is  also  a 
lack  of  sensitivity  to  elevation  angles  near  broadside  for  a  2-element  antenna  array.  To 
accurately  determine  the  angle  of  arrival  of  a  wave  both  in  azimuth  and  elevation,  the 
phase  measurements  of  two  or  more  interferometer  baselines  must  be  combined. 
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2.  Multiple  Baselines 


Interferometers  with  more  than  one  baseline  allow  the  use  of  longer  baselines 
to  increase  the  aperture  of  the  array  while  keeping  at  least  one  short  baseline  (less  than 
one-half  wavelength)  to  resolve  ambiguities. 

By  increasing  the  aperture  of  the  array,  the  angle  of  arrival  of  a  signal  can  be 
determined  with  greater  precision.  The  resolution  of  a  multiple  baseline  interferometer  is 
given  by  the  widest  antenna  pair  (baseline). 

This  study  will  analyze  a  multiple  baseline  phase-interferometer  antenna  array 
developed  by  SWRI.  In  this  interferometer  the  3-dimensional  angle  of  arrival  of  a  plane 
wave  is  resolved  by  using  the  phase  across  two  apertures  that  are  not  parallel. 

Consider  the  coordinate  system  for  a  non-parallel  baseline  interferometer,  as 
shown  in  Fig.  4,  with  three  antennas  located  at  P^(0,0,0),  F,(A:,,y,,2,),  and  P^{x^,y^,z^) . 

Antenna  Po  is  the  reference  antenna  and  its  phase  center  is  located  at  the 
earth’s  surface.  The  computed  values  of  azimuth  and  elevation  will  be  referenced  to  a 
plane  which  is  tangent  to  the  earth’s  surface  at 

Let  the  shortest  distance  between  a  normal  to  the  plane  wave  that  goes  through 
antenna  and  the  point  P,  be  Similarly,  let  the  shortest  distance  between  a 

nonnal  to  the  plane  wave,  passing  through  antenna  P^  and  the  point  P^  be  P^, . 
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Fig.  4.-  Coordinate  System  for  Non-Parailei  Baseline  Interferometer. 


Assuming  a  plane  wave  arrivitig  from  azimuth,  (|> ,  and  elevation,  y ,  where 
y  =  90°  -  0,  the  phase  difference  induced  on  each  base  line  is  given  by 

Phase  difference  PJP^  -  P(?^)  ®  ^ot 


Phase  difference  P/j  = 

where 

PJ^j^  =  distance  from  P^  to  and 
PJP^  -  distance  from  P^  to  P^ . 

To  derive  an  expression  for  the  coordinates  of  point  P^  must  first  be 
found.  From  Fig.  4, 
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(10) 


therefore 


in  addition 


then 


and  finally 


The  coordinates  of  can  now  be  expressed  as 


Xj^  tan(|)  , 


cos<]> 


(11) 


(12) 


(13) 


(14) 
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and 


Next,  to  minimize  this  distance 


/ 


- -2. 

cos<|> 


dCFpy 

=  0  . 
dx 

Using  this  expression  and  solving  for  and  yields 


Xj^  =  jc,  cos*^  cosV‘^)’i  cos*Y+2,  cos(l>  cosYsiny  , 


and 


2^  *  x^sec(|)tanY  . 


(17) 


(18) 


(19) 


(20) 


(21) 


The  distance  can  then  be  found  by 


'’A  ” 

(22) 

or 

P^^  =  jc,  cos<|»  COSY  *  y,  sin<])cosY  +  2,  sinY  . 

(23) 

The  same  procedure  can  be  repeated  for  PJP^  to  find 
PJP^  =  JCj  cos(j)  COSY  *  >’2  h 
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For  a  plane  wave  arriving  from  azimuth,  <|> ,  and  elevation,  y,  it  is  possible 
to  express  the  phase  difference  induced  in  each  baseline  measurement  as  the  distance  a 
wave  travels  between  points  and  to  point  as 

P  =  __  (x,  cos(|»  cos\j/  +  y^  sin()>  cosy  +  r,  sin\|r)  (25) 

\ 

and 

Pqj  =  —  (JCj  cos(j)  cos\|/  +  sin(j>  cos\|/  +  Z^sinv;/)  .  (26) 

X 

In  practice,  the  values  of  the  coordinates  of  tlie  three  antennas  and  the 
frequency  of  the  incoming  signal  are  known,  and  the  phase  difference  P^i  ^02 

be  measured.  The  two  equations  and  two  unknowns  can  be  solved  for  the  azimuth,  <}> ,  and 
elevation,  y,  of  the  incoming  wave. 

To  simpliiy  the  process,  the  following  variables  can  be  defined  as 
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and 


^  . 

c  =  y,x^  -  . 


(29) 

(30) 

(31) 

(32) 


C.  IONOSPHERIC  CONSIDERATIONS 
1.  Propagation  Modes 

The  ionosphere  plays  an  essential  role  in  the  determination  of  target  location 
when  using  an  HFDF  SSL  system.  The  ionosphere  is  the  region  of  the  atmosphere 
composed  of  gases  which  are  primarily  ionized  by  solar  radiation  in  the  ultra-violet  to 
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X-ray  range.  It  extends  from  about  50  km  to  1000  km  in  altitude  and  is  an  amorphous, 
nonhomogeneous,  and  extremely  turbulent  plasma  [Ref.  3]. 

The  ionosphere  is  usually  modeled  as  a  series  of  layers  that  vary  in  height  and 
electron  density.  For  all  practical  purposes  3  layers  have  been  identified  in  the  ionosphere 
D,  E,  and  F. 

The  D-layer  extends  from  50  km  to  90  km  and  exists  only  during  daylight 
hours.  It  is  an  absorptive  layer,  where  the  lower  HF  frequencies  are  attenuated  more  than 
higher  frequencies. 

The  E-layer  extends  from  90  km  to  130  km  and  although  at  night  it  is  quite 
weak,  it  is  considered  a  permanent  layer.  During  daylight  it  is  used  for  short-to-medium 
range  propagation. 

In  the  E-region  height  at  about  100  km,  a  highly  ionized  layer  occurs  at 
various  times  known  as  the  Sporadic  E  (Es)  layer.  This  transient  layer  of  short  duration 
appears  during  daylight  allowing  unexpected  propagation  of  high  frequencies. 

The  F-layer  is  the  most  important  layer  for  long  distance  communications.  It 
exists  at  130  and  500  km  altitude  and  is  a  permanent  layer.  Throughout  the  year  in  low 
latitudes  and  during  the  summer  in  high  latitudes,  the  F  region  splits  into  two  layers,  F, 
and  Fj,  after  sunrise. 

The  lower  layer,  F„  is  at  a  height  of  approximately  200  km  and  is  present  only 
during  daylight  hours.  The  higher  layer,  Fj,  occurs  at  approximately  300  km  and  has  the 
highest  ionization  density. 

A  plot  of  ionospheric  electron  density  versus  height  is  shown  in  Fig.  5. 
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Fig.  5.-  Electron  Density  of  the  Ionosphere. 

During  daylight,  the  lower  frequencies  of  the  HP  band  are  highly  attenuated 
by  the  D-layer  and  efficient  long  distance  communication  is  in  the  range  of  10  to  20 
MHz.  At  night,  higher  frequencies  are  not  reflected  from  the  weaker  nighttime  ionosphere 
and  the  lower  portion  of  the  HP  band  (2  to  10  MHz)  must  be  used. 

To  comply  with  the  tactical  scenario  the  HPDP  SSL  system  must  then  be  able 
to  cover  2  to  20  MHz. 

The  term  "propagation  mode"  refers  to  the  path  of  the  signal  through  the 
ionosphere  and  the  number  of  reflections  (hops)  that  occur  between  the  transmitter  and 
receiver  (emitter  to  HPDP  SSL  in  this  case).  Propagation  modes  are  identified  by 
numbers  and  letters  which  represent  die  number  of  hops  and  the  name  of  the  layer  where 
each  reflection  occurs. 
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Examples 

2F2  =  A  path  with  2  reflections  from  the  Fj-layer. 

FjEs  =  One  reflection  from  the  Fj-layer  followed  by  a  reflection  from  the 
Sporadic  E-layer  (also  called  an  N-type  reflection). 

FjEsFj  =  One  reflection  from  the  Fj-layer  followed  by  a  reflection  from  the 

Sporadic  E-layer  and  then  another  Fj-hop  (called  an  M-type  reflection). 

Proper  identification  of  the  mode  of  propagation  is  critical  for  HFDF  SSL 
operation.  Incorrect  path  assumptions  create  errors  in  final  position  determination. 

2.  Ionospheric  Tilts  and  other  Effects 

Variations  of  the  ionization  gradient  with  geographical  position  are  usually 
called  "ionospheric  tilts".  Because  the  elective  tilt  of  an  ionospheric  layer  depends  on  the 
depth  of  penetration  of  a  ray  into  the  layer,  a  tilted  layer  should  not  be  considered  as  a 
tilted  mirror  [Ref.  2].  The  tilt  depends  on  the  layer  as  well  as  the  ray  path  and  varies  from 
ray  to  ray.  Ionospheric  tilts  can  be  lateral  or  longitudinal  with  respect  to  the  path  between 
the  emitter  and  the  HFDF  SSL. 

Several  methods  are  proposed  for  obtaining  information  on  ionospheric  tilts 
and  for  correcting  their  effects  on  azimuth  and  elevation  angles  measured  at  an  HFDF 
SSL.  The  Chapman  model  and  numerical  maps  based  on  monthly  median  predictions  of 
ionospheric  parameters  are  two  methods  that  do  not  require  real-time  measurements.  The 
use  of  a  "check"  transmitter  of  known  position,  the  direct  measurement  of  local  tilt  by  a 
pulsed  sounder  and  doppler  measurements  or  ionospheric  profiles  from  vertical  or  oblique 
incidence  sounders  provide  real-time  information  on  tilts.  A  more  detailed  treatment  of 
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each  of  these  methods  is  found  in  Ref.  2.  The  DF  system  in  this  study  uses  a  vertical 
incidence  sounder  to  determine  tUt  characteristics  of  the  layer  of  interest. 

For  a  single  isolated  propagation  mode,  the  ionosphere  acts  as  a  mirror  and  the 
spread  of  rays  within  a  given  mode  is  quite  small.  When  complex  mode  structures  are 
present  they  produce  large  number  of  ray  paths.  The  received  signal  is  then  the  vector 
sum  of  randomly  varying  ray  components  and  the  resulting  signal  intensity  varies 
randomly.  This  effect  is  known  as  fading  and  is  present  in  both  short  and  long  range 
propagation  paths,  but  for  this  study  only  the  300  to  500  km  paths  are  considered.  In 
addition  to  fading,  sky  waves  can  exhibit  changes  in  direction  of  arrival,  both  in  azimuth 
and  in  elevation  [Ref.  4].  These  changes  can  be  rapid  fluctuations  generally  attributed  to 
lack  of  homogeneity  of  the  ionosphere,  or  long-term  fluctuations  due  to  large  scale 
ionospheric  disturbances.  Small  variations  in  the  characteristics  of  the  ionosphere  can  arise 
from  traveling  ionospheric  disturbances  (TID)  and  ISrom  other  temporal  irregularities. 

Time  averaging  can  be  used  to  counteract  signal  variations  and  to  resolve  the 
dominant  propagation  mode  when  using  a  phase-interferometer  HFDF  SSL.  The  DF 
system  of  this  study  uses  a  sliding  window  of  20  frames  which  is  the  average  observation 
over  a  one-second  period  [Ref.  5]. 

The  earth’s  magnetic  field  is  another  important  factor  to  consider  in  HFDF 
SSL  operation  since  it  modifies  the  phenomena  of  ionospheric  refraction.  It  "splits"  an 
incident  wave  into  two  polarization  components,  each  with  different  propagation  velocities 
and  attenuations.  These  components  are  called  the  ordinary  and  extraordinary  waves.  The 
extraordinary  wave  suffers  much  higher  absorption,  thus  useful  sky  wave  propagation 
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almost  always  takes  place  through  the  ordinary  wave  pief.  6].  For  transverse  or  quasi- 
transverse  propagation,  the  usual  mode  for  HF  sky-wave  propagation  except  near  the 
magnetic  poles,  the  critical  frequency  of  the  ordinary  wave  is  the  same  as  that  calculated 
neglecting  the  earth’s  magnetic  field  [Ref.  6].  Only  the  reflected  ordinary  wave  must  be 
considered  in  this  study. 

The  earth’s  magnetic  field  also  changes  the  polarization  of  the  reflected  wave 
with  respect  to  the  incident  wave.  A  plane-polarized  incident  wave  is  converted  into  an 
elliptically-polarized  reflected  wave.  When  the  direction  of  propagation  forms  an  acute 
angle  with  the  direction  of  the  magnetic  field,  the  reflected  ordinary  wave  will  be  left- 
hand  elliptically  polarized,  and  the  extraordinary'  wave  right-hand  elliptically  polarized. 
If  the  angle  is  obtuse,  the  reverse  case  is  true.  Generally  in  the  northern  hemisphere  the 
downcoming  ordinary  wave  is  polarized  in  the  left-hand  sense  and  in  the  southern 
hemisphere  it  is  polarized  in  the  right-  hand  sense  [Ref.  6].  The  incoming  sky  wave  will 
then  have  vertical  and  horizontal  electrical  field  components  and  since  the  interferometer 
compares  the  signal  received  by  two  identical  antennas,  the  polarization  of  the  antennas 
is  not  a  critical  issue.  Anteima  polarization  can  become  an  important  issue  if  both  receive 
antennas  are  not  correctly  aligned  with  each  other  or  if  the  polarization  characteristics  are 
not  the  same. 

An  extensive  analysis  of  polarization  error  of  a  two  element  interferometer  has 
appeared  in  the  Russian  literature  [Ref.  7].  These  studies  show  that  significant  phase 
errors,  approaching  180°,  can  occur  independent  of  spacing  error,  when  the  antennas  do 
not  have  identical  polarization  characteristics. 
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D.  HFDF  SSL  VERTICAL  TRIANGULATION 


The  goal  of  an  HFDF  SSL  is  to  detennine  the  position  of  a  transmitter  with  a 
moderate  degree  of  accuracy  and  over  a  reasonably  wide  range  of  distances,  operating 
frequencies  and  ionospheric  conditions.  From  measurements  of  die  angle  of  arrival  of  a 
signal  position  ^ing  is  achieved  using  a  technique  called  vertical  triangulation.  Utilizing 
a  single  station  makes  the  system  simple  and  inexpensive  to  install  and  maintain, 
compared  to  a  network  of  DF  stations. 

Appleton  and  Barnett  first  used  vertical  triangulation  in  HF  propagation  in  1924 
[Ref.  8],  using  the  geometry  of  Fig.  6.  They  applied  a  thin-ionosphere  approximation 
model,  assuming  a  mirror  reflection  at  point  L.  From  the  distance,  d,  and  elevation  angle 
of  arrival,  y,  they  deduced  the  height  of  the  ionosphere. 

This  procedure  also  assumes  a  uniform  ionosphere,  a  fiat  earth  and  straight  ray 
paths.  Later  studies  have  demonstrated  that  the  structure  of  the  ionosphere;  is  not  uniform. 
Refraction  bends  ray  paths  and  the  efiect  of  curvature  of  the  earth  must  be  included  for 
long  distances. 

The  HFDF  SSL  problem  of  emitter  location  is  solved  from  receiver  site  information. 
Knowing  the  elevation  and  azimuth  angles  of  arrival  of  the  signal  and  the  height  of  the 
ionosphere,  the  bearing  and  distance  to  the  emitter  are  determined.  For  accurate 
measurements,  the  effects  of  ray  refraction  in  the  atmosphere,  earth  curvature  and 
ionospheric  structures  should  be  included.  Ray  path  determination  is  simple  up  to  the 
point  where  rays  enter  the  ionosphere.  The  exit  path  from  the  ionosphere  to  the  HFDF 
SSL  station  depends  on  the  electron  density  distribution. 
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IONOSPHERE 


h - H 

d 

Fig.  6.-  Vertical  Triangulation  Using  Appleton’s  Model. 

If  this  distribution  were  accurately  known,  ray  tracing  methods  for  calculating 
emitter-receiver  distance  are  available  a?  computer  codes.  Unfortunately,  the  electron 
density  structure  of  the  ionosphere  is  rarely  known  with  sufficient  precision  to  determine 
the  amount  of  ionospheric  penetration  which  has  occurred. 

A  vertical  sounder  can  be  used  to  determine  the  virtual  height  (h’)  of  the  apparent 
reflection  point  as  shown  in  Fig.  7.  With  the  virtual  height  and  the  angle  of  arrival  of  the 
incoming  signal,  vertical  triangulation  is  used  to  calculate  distance  to  the  emitter. 

The  HFDF  SSL  system  under  study  uses  a  sounder  to  obtain  ionospheric 
information  for  the  tactical  scenario  developed  in  Chapter  I.  With  300  to  500  km  ranges 
expected,  vertical  triangulation  can  be  performed  using  the  flat-earth,  flat-ionosphere 
approximation  [Ref.  9]. 
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Fig.  7.-  Signal  Geometr}’  for  the  Flat  Earth  Flat  Ionosphere  Approximation. 


The  distance  to  the  emitter  is: 


d  =  ,  (35) 

tan  \|/ 

where 

h’  =  virtual  height  of  the  ionosphere, 
d  =  distance  between  emitter  and  receiver  and 
\}r  =  elevation  angle  of  arrival  of  the  signal. 

From  this  expression  a  family  of  curves  of  virtual  height  versus  elevation  angle  for 
various  distances  are  produced  in  Fig.  8. 

E.  TYPICAL  HFDF  SSL  PATHS 

Based  on  ionospheric  considerations  of  this  chapter  and  the  tactical  scenario  of 
Chapter  I,  the  characteristics  of  typical  HFDF  SSL  paths  can  be  established. 
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Fig.  8.-  Transmission  Curves  Parametric  in  Distance  (d)  for  Flat-Earth  Flat-Ionosphere. 


The  objective  of  this  research  is  to  analyze  the  effects  of  errors  introduced  in  the 
HFDF  SSL  process  by  the  antenna  array.  For  this  investigation  the  assumptions  are  (1) 
the  virtual  height  of  the  ionosphere  is  known  and  (2)  the  effects  of  ionospheric  tilts  can 
be  accounted  for. 

To  determine  the  performance  of  the  antenna  arrays,  values  of  the  minimum  and 
maximum  elevation  angles  of  arrival  of  the  incoming  waves  must  be  determined.  Using 
the  transmission  curves  of  Fig.  8,  for  virtual  heights  of  the  E-layer  and  the  Fj-Iayer  at 
distances  between  300  and  500  km,  it  is  seen  that  elevation  angles  of  anival  vary  from 
25°  to  65°. 
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in.  NEC  MODELING 


A.  BACKGROUND 

The  Numerical  Electromagnetics  Code  (NEC)  is  designed  for  modeling  of  the 
electromagnetic  response  of  general  structures.  NEC  was  developed  by  the  Lawrence 
Livermore  National  Laboratory  for  the  Naval  Ocean  Systems  Center  (NOSC).  It  uses  an 
electric  field  integral  equation  to  model  wire-like  objects,  and  a  magnetic  field  integral 
equation  to  model  closed  surfaces  with  time  harmonic  excitation  [Ref.  10].  For  stmctures 
containing  both  wires  and  surfaces  the  two  equations  are  coupled  together. 

NEC  solves  the  integral  equations  numerically  via  the  Method  of  Moments.  Finitely 
conducting  ground  effects  can  be  included  in  the  model,  in  addition  to  free  space  and 
perfect  ground  conditions.  Two  alternative  models  are  provided  for  finitely  conducting 
ground.  The  Reflection  Coefficient  Approximation  uses  image  fields  modified  by  Fresnel 
reflection  coefficients  and  is  valid  for  structures  located  at  least  0.1  wavelengths  above 
ground.  The  Sommerfeld  integral  formulation  for  the  near  field  at  the  interface  provides 
an  accurate  solution  for  elements  above,  below,  or  penetrating  the  ground,  but  with  a 
substantial  cost  in  computational  resources. 

NEC  permits  excitation  by  voltage  sources  on  the  structure  or  by  incident  plane 
waves  of  specified  polarization.  Outputs  include  radiation  patterns,  power  and  directive 
gains,  feed  point  impedances,  ohmic  losses,  efficiency,  wire  and  surface  currents  and 
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charge  density,  maximum  couplings,  near  E  and  H  fields,  scattering  cross  section,  and 
average  gains.  A  detailed  description  of  the  code  is  found  in  Ref.  11. 

B.  SWRI  HFDF  SSL  SYSTEM  DESCRIPTION 

The  HFDF  SSL  system  modeled  in  this  thesis  employs  the  phase-interferometer 
technique.  A  two  orthogonal-baseline  antetma  array  determines  the  three-dimensional 
angle  of  arrival  ol  a  signal  as  described  in  Chapter  n.  The  system  has  been  developed  by 
SWRI  [Ref.  12],  and  is  designed  to  be  a  real-time  computer-controlled  HFDF  SSL.  It 
provides  DF  capabilities  for  ground  and  sky  wave  HF  propagation  for  the  2  to  30  MHz 
frequency  range  at  ranges  up  to  1 100  km.  The  system  is  intended  for  tactical/transportable 
operation.  The  system  block  diagram  is  shown  in  Fig.  9.  Three  main  components  can  be 
identified  as; 

(1)  two,  seven-element  interferometer  antenna  arrays, 

(2)  automatic  DF  acquisition,  processing  and  display  system  and 

(3)  ionospheric  sounder  system. 

This  research  study  will  concentrate  on  the  analysis  of  the  two  antenna  arrays.  The 
DF  processing  algorithm  used  by  SWRI  is  used  in  the  evaluation  of  the  effects  of  errors 
induced  by  the  antenna  arrays  under  different  ground  conditions. 

1.  Low  Band  Antenna  Array 

The  low  band  antenna  array  is  intended  to  cover  2  to  10  MHz  by  seven 
identical  vertical,  square  crossed  loops.  The  antenna  elements  are  placed  in  an  "L"  shaped 
array  as  shown  in  Fig.  10.  Each  orthogonal  aim  of  the  array  contains  three  antennas 
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which  are  placed  along  a  line  at  different  distances.  A  central  antenna  (reference)  is 
located  on  the  vertex  of  the  array  and  allows  multiple  baselines  to  be  selected  on  the 
interferometer.  It  is  diagonally  symmetric,  that  is,  each  baseline  on  one  atm  (formed  by 
the  central  reference  antenna  and  one  of  the  antennas  of  the  atm)  has  an  image  on  the 
other  arm.  Because  of  this  characteristic  the  system  can  compare  the  phase  of  one 
baseline  with  its  equivalent  on  the  other  arm  and  by  applying  phase-interferometer 
techniques  it  can  find  the  angle  of  arrival  of  a  signal. 

The  shortest  baseline  of  the  low  band  array  is  13  meters  long,  less  than  one- 
half  wave  length  at  10  MHz  (15  meters).  This  is  the  highest  ftequency  of  the  array  for 
unambiguous  operation.  The  characteristics  of  the  low  band  elements  are: 

Antenna  type:  Crossed  Loops 

Orientation:  Vertical 

Element  shape:  Squaie  loop 

Frequency  range:  1.5  -  10  MHz 

Size:  5  foot  /  side 

Number  of  elements:  2  orthogonal  loops 

Height  above:  0.15  m 

Vertical  distance  between  elements:  0.05  m 

Cross  section:  0.05  m 

Material:  Aluminum  and  fiberglass 


The  geometry  of  a  single  element  of  the  array  is  detailed  in  Fig.  11. 


2.  High  Band  Antenna  Array 

The  high  band  antenna  anay,  designed  for  10  to  30  MHz  has  seven  identical, 
vertical  monopoles,  placed  over  a  wire  mesh  ground  screen  (counterpoise).  It  is  also  an 
"L"  shaped  array  but  with  different  dimensions  than  the  low  band  one.  It  is  also 
diagonally  symmetric  as  shown  in  Figure  12. 
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Fig.  11-  5  Ft  Square  Crossed  Loop. 


_ 50  m 

33.33  re 


(Drawing  not  to  scale) 


Fig.  12.-  Interferometer  High  Band  Antenna  Ajrray. 


The  shortest  baseline  of  the  array  is  4.33  m,  and  is  used  to  resolve  ambiguities. 
This  is  less  than  one-half  wavelength  at  30  MHz,  the  highest  frequency  of  operation.  The 
characteristics  of  the  high  band  array  elements  are 
Antenna  type:  Monopole 

Orientation:  Vertical 

Frequency  range:  10  -  30  MHz 

Size:  6  feet 
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Height  over  ground: 


0-15  m 


Cross  section: 
Material: 

Ground  screen: 
Screen  radials: 

Radial  wire  diameter: 


0.05  m 

Aluminum  fiberglass 
6x6  feet  grid  with6  inches  spacing 
6  feet  radials  every  45® 

0.01  m 


The  shape  of  a  single  element  of  the  high  band  antenna  array  is  shown  in 

Fig.  13. 

C.  COMPUTER  MODEL  DEVELOPMENT 
1.  Dimensions  Considerations 

The  first  step  in  developing  a  numerical  model  of  the  two  antenna  arrays  is  to 
define  the  geometry  of  the  elements  and  their  physical  location  within  the  array. 

The  second  step  is  selecting  the  ’mgth  of  wire  segments  used  to  model  each 
antenna  structure  within  the  guidelines  of  the  NEC  code  specifications.  Since  cunent  is 
computed  at  the  center  of  each  segment,  the  length  of  the  segments  partially  determines 
the  resolution  of  the  currents  in  a  NEC  model.  Critical  regions  of  the  antenna  require 
smaller  segment  lengths. 

The  maximum  segment  length  recommended  in  NEC  is  <0.1  A .  Additionally, 
the  NEC  segment  length  for  32-bit  double  precision  applications  is  >  10  *  A . 
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Fig.  13."  6  Ft  Monopole  With  6  Ft  Sq.  Wire  Mesh  and  Radials. 
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Applying  these  two  constraints  to  2-20  MHz  the  segment  length  limits  were 
calculated  as: 

Frequency  Wavelength  Maximum  Minimum 


2 

MHz 

150 

m 

< 

15  m 

0.15 

A 

10"’ 

m 

5 

MHz 

60 

m 

< 

6  m 

0.6 

A 

10"® 

m 

10 

MHz 

30 

m 

< 

3  m 

0.3 

A 

10"® 

m 

15 

MHz 

20 

m 

< 

2  m 

0.2 

A 

10"® 

m 

20 

MHz 

15 

m 

< 

1.5  m 

0.15 

A 

10"® 

m 

From  these  calculations  for  the  worst  case  conditions,  the  segment  lengths 
should  be  within  the  range  of  0.15  x  10*’  m  to  3  m  for  the  low  band  array  and  0.3  x  lO"* 
m  to  1.5  m  for  the  high  band  array.  The  model  was  adjusted  to  verify  that  it  contained 
the  minimum  number  of  segments  to  reduce  computation  time,  whUe  maintaining  the 
segment  lengths  within  the  limits  of  acceptable  calculation  accuracy.  The  segment  lengths 
of  the  model  were  selected  to  be  from  0.51  m  to  1.72  m  for  low  band  and  from  0.15  m 
to  0.91  m  for  high  band.  In  addition,  segment  lengths  were  selected  to  comply  with  a 
NEC  restriction  for  errors  of  less  than  1%  as 

Segment  length  /  wire  radius  >  2. 

Segment  length  limits  had  previously  been  selected  so  dimensions  of  the  radius 
of  all  the  wires  in  the  model  were  held  to 

low  band  <  0.25  m  and 
high  band  <  0.07  m. 


39 


The  maximum  wire  radius  used  in  the  computer  model  was  0.05  m  for  both 
arrays,  less  than  the  NEC  requirement  above.  Also,  when  modeling  structures  over  ground 
NEC  imposes  a  restriction  for  wires  placed  horizontally  close  to  ground  as 

*  a^)^>  10-«A.  ,  (36) 

where 

h  =  height  of  the  horizontal  wire  over  ground, 
a  =  wire  radius  and 
X  =  wavelength. 

Only  the  high  band  antetuia  array  contained  horizontal  wires  close  to  ground 
(the  wire  mesh).  For  these  conditions  and  20  MHz,  the  worse-case  frequency,  the  wire 
height  of  the  mesh  was  set  0.001  m  above  ground  to  satisfy  NEC  guidelines  for  horizontal 
wires  near  ground. 

2.  Symmetry  Considerations 

NEC  includes  capability  for  exploitation  of  structure  symmetry  to  simplify 
model  development  and  to  reduce  computational  requirements.  Both  arrays  exhibit  partial 
diagonal  symmetry  and  by  using  the  Numerical  Green’s  Function  (NGF)  option  of  NEC, 
the  symmetric  portions  of  the  stractures  can  be  evaluated  and  the  factored  interaction 
matrix  saved  in  a  file  for  later  use.  The  complete  solution  can  be  obtained  without  repeat¬ 
ing  calculations  for  the  symmetric  data  on  file.  This  feature  is  particularly  useful  for  the 
high  band  model  because  of  the  large  number  of  segments  involved  in  modeling  the  wire 
mesh  for  every  monopole  of  the  array. 
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3.  Ground  considerations 


The  Sonunerfeld  method  of  ground  interaction  was  used  to  provide  an  accurate 
solution  for  the  wires  close  to  ground  in  both  arrays.  "Poor"  ground  and  "good"  ground 
parameters  were  used  in  addition  to  free  space  and  perfect  ground  to  thoroughly  cover 
effects  which  were  influenced  by  installation  conditions.  The  ground  constants  chosen  are 


the  following 

GROUND 

CONDUCTIVITY 

RELATIVE  DIELECTRIC 

CONDITION 

(O) 

CONSTANT  (E^) 

Good  Ground 

0.01  mhos/m 

30 

Poor  Ground 

0.001  mhos/m 

5 

These  ground  parameters  were  used  with  a  separate  program  called  SOMNTX 
to  generate  a  ground  interaction  data  file,  containing  values  required  by  NEC  when  using 
the  Sommer-feld  option. 

D.  NEC  DATASETS 

The  NEC  datasets  constructed  for  the  high  and  low  band  antenna  arrays  used  by  the 
HFDF  SSL  phase-interferometer  are  listed  in  Appendix  A. 
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IV.  PERFORMANCE  OF  THE  ANTENNAS 


An  HFDF  SSL  system  is  always  used  to  analyze  the  psrformance  of  the  antennas 
in  the  receive  mode  so  the  Reciprocity  Theorem  was  used  to  characterize  the  arrays  for 
all  performance  parameters  except  feed  point  receive-current  phases.  When  calculating 
angle  of  arrival,  the  more  resource-consuming  receive  mode  of  excitation  had  to  be  used. 

A.  ANTENNA  INPUT  IMPEDANCE 

Input  impedance  is  an  important  performance  parameter  of  an  antenna  and  because 
its  value  depends  on  the  electrical  configuration  and  dimensions  of  the  antenna  it  provides 
an  indication  of  the  validity  of  a  numerical  antenna  model. 

To  verify  the  accuracy  of  the  crossed-loop  NEC  model  the  impedance  from  NEC 
was  compared  to  values  calculated  by  a  theoretical  method. 

In  the  2  to  10  MHz  frequency  range,  the  5  foot  square  loops  are  a  small  antenna 
compared  to  the  wavelength.  These  small  loops  behave  as  if  they  were  "magnetic  dipoles" 
with  liighly  inductive  input  impedance.  The  method  proposed  by  Ref.  13  was  used  to 
calculate  the  input  impedance  of  the  vertical  square  loops.  This  method  assumes  that  the 
square  loop  can  be  analyzed  as  a  shorted  transmission  line  of  wire  radius  which  is  very 
small  compared  to  the  loop  dimensions  and  to  the  wavelength.  The  reactive  part  of  the 
input  impedance  is  a  function  of  the  loop  perimeter  wire  radius.  Thus  an  equivalent 
transmission  line  shorted  at  one  end  can  be  used  to  determine  its  impedance.  Using  this 
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model,  the  input  impedance  of  the  shorted  lossless  transmission  line  is  equal  to  the 
reactive  part  (Xin)  of  the  input  impedance  of  the  loop  as 

X,^^jZ„tm{kl)  ,  (36) 

where 

k  =  wave  number  or  propagation  coefficient  and 
1  =  transmission  line  length. 


Without  losing  generality,  the  transmission  line  length,  1,  can  be  approximated  as 
half  the  perimeter  and  the  characteristic  impedance  of  the  transmission  line  (Z^)  then 


becomes 


Z,  =  276  log 


s 


(37) 


where 

1  =  (loop  perimeter )/2, 
s  =  (loop  area)/l  and 
r  =  wire  radius  of  the  loop. 

Applying  the  preceding  fonnulas  at  5  MHz,  a  wire  radius  of  0.05  m,  and  a  loop  side 
length  of  1 .524  m. 


/ 


{  1.524  X  4 


■>1 

/ 


3.048  m  , 


(38) 


s 


1.524"  "I 
3.048  J 


0.762  m 


(39) 
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and 


k 


^  2n  ^ 

"  2jt  ' 

.  ^  J 

[60  J 

0.104  m  . 


Then 


Z.  =  276  log 


0.76 

0.05 


=  326.5  Ohms 


(40) 


(41) 


and 


X.„  =  j  (  326.5  )  tan  [  0.104  x  3.04  ] 
=  j  107.9  Ohms. 


(42) 


The  NEC  value  for  the  imaginary  part  of  the  input  impedance  was  110.5  Ohms, 
which  agrees  well  with  the  approximate  theory. 

For  the  resistive  part  of  the  loop  input  impedance  the  least  squares  technique  of 
Ref.  13  yields 


=  a 

tan 

'  kp 

V 

2 

• 

K 

(43) 


where 

k  =  wave  number  or  propagation  coefficient  and 
p  =  perimeter  of  the  loop. 
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In  addition,  a  =  1.126  and  b  =  3.95  for  a  square  loop.  Using  this  formula,  Ri„  = 
0.014  Ohms.  The  resistive  part  of  the  input  impedance  from  a  NEC  calculation  agrees 
well  at  0.013  Ohms. 

The  high  band  antenna  array  is  composed  of  seven  vertical  monopoles,  each  one 
with  an  underlying  wire  mesh.  To  verify  the  accuracy  of  the  NEC  model,  the  input 
impedance  of  a  monopole  of  the  same  characteristics  was  calculated  and  compared  to  the 
NEC  result. 

The  input  resistance  of  an  ideal  monopole  over  perfect  ground  is  half  the  radiation 
resistance  of  a  dipole  whose  length  is  twice  the  monopole  height.  From  Ref.  14  the  radia¬ 
tion  resistance  of  a  short  dipole  is 

IP  (  T  w  ^ 

R  ^  =  20  _  .  (44) 

I  /„  P  I  j 

Since  the  height  of  the  vertical  monopole  is  1.83  m,  the  equivalent  length  of  its 
dipole  counterpart  is  3.66  m.  By  using  this  length  (L)  and  a  frequency  of  15  MHz, 

R  .  =  20  f  T  =  6.6  Ohms  .  (45) 

20 

The  radiation  resistance  of  the  monopole  is  thus  approximately  3.3  Ohms,  assuming 
a  lossless  case.  The  value  of  the  input  resistance  of  the  monopole  calculated  by  NEC  code 
was  3.74  Ohms,  close  enough  to  be  considered  adequate. 
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B.  SMITH  CHART  INPUT  IMPEDANCE  ANALYSIS 


An  unportant  verification  of  the  behavior  of  an  antenna  through  a  frequency  range, 
is  to  analyze  its  input  impedance  using  the  Smith  Chart.  When  plotted  on  a  Smith  Chart, 
the  input  impedance  of  the  antennas  should  show  smooth  variation  with  clockwise 
rotations  and  describe  a  smooth  curve,  as  frequency  is  increased  over  the  operating  range. 

The  model  of  the  electrically  small  crossed  loop  was  first  tested  for  the  2  to  5  MHz 
range  for  good  ground  (e=30  and  o=0.01),  poor  ground  (e=5  and  a=0.001),  perfect 
ground  and  free  space  conditions.  The  results  of  the  different  runs  are  shown  in  Table  B-1 
of  Appendix  B.  Each  set  of  data  was  then  nomialized  to  the  lowest  value  of  its  set  and 
the  resultant  input  impedances  plotted  on  a  Smith  Chart. 

As  shown  in  Fig.  14  and  Fig.  15,  the  input  impedfuices  show  smooth  variation  and 
describe  highly  reactive  curves  on  the  Smitli  Chart,  confirming  that  the  NEC  results  of 
the  small  crossed  loop  impedance  are  as  expected. 

Next,  the  NEC  model  of  the  vertical  monopole  with  wire  mesh  was  used  to 
construct  a  Smith  Chart  of  input  impedance  for  10  to  20  MHz  for  free  space,  good 
ground,  poor  ground,  and  perfect  ground  conditions.  The  results  are  in  Table  B-2  of 
Appendix  B.  The  resulting  Smith  chart  plots  are  shown  in  Fig.  16  and  Fig.  17  and  are  as 
expected. 

C.  AVERAGE  POWER  GAIN 

Another  parameter  used  to  verify  numerical  antenna  models  is  the  average  power 
gain  of  the  antenna.  Since  it  accounts  for  antenna  efficiency,  hence  radiation  resistance, 
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Fig.  17.*  Monopole  Input  Impedance  (Good  and  Poor  Ground). 


it  is  a  measure  of  validity  of  current  distribution  at  the  feed  point. 


The  average  power  gain  is  defined  as  the  integral  of  the  antenna  power  gain  per  unit 
solid  angle, 

-  z  /o 

where 

^avg  ~  average  power  gain, 

=  antenna  power  gain  and 
Q  =  solid  angle. 

When  the  solid  angle  is  4 n,  the  average  gain  can  be  expressed  as 


where  is  the  total  power  radiated,  Pi„  is  the  input  power  at  the  antenna  terminals, 


and  K  is  1  for  free  space,  2  for  perfect  ground  and  between  0  and  1  for  lossy  ground. 


can,  in  turn,  be  found  by  integrating  the  complex  pointing  vector,  5 


as 

jj{ExH-)ds  ,  (48) 

where 

E  =  electric  field  intensity  and 
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H'  =  complex  conjugate  of  the  magnetic  field  intensity. 


In  the  far  field, 

P  =  .£!  r  J_£_L  d  Cl  for  R  CO  ,  (49) 

2  r\ 

where 

T|  =  intrinsic  impedance  of  the  medium  and 
R  =  distance. 

The  input  power  at  the  antenna  terminals  can  be  found  by 

/>,.  =  i  SeinnO)]  ,  (50) 


where 

Kj,  =  source  voltage  and 

/*(0)  =  complex  conjugate  of  current  at  antenna  temunals. 

For  a  lossless  antenna  the  radiated  power  is  equal  to  the  input  power  and  the 
efficiency  is  100%.  The  average  gain  will  depend  upon  K  (i.e.  the  ground  conditions). 
The  lossless  HFDF  SSL  antenna  elements  were  first  considered  in  free  space  where  the 
expected  value  of  the  average  power  gain  is  1.  For  operation  over  ground,  average  power 
gain  should  be  2.  The  NEC  model  was  also  tested  over  the  two  types  of  lossy  ground 
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specified  previously,  good  ground,  £=30  and  a=0.01,  and  poor  ground,  £=5  and 


C=0.001. 

The  average  gains  for  the  NEC  models  of  the  crossed  loop  and  monopole  antennas 
with  these  ground  conditions  are  shown  in  Figures  18  and  19 

The  average  gains  for  the  NEC  models  are  well  within  ±  10%  of  theoretical  values 
for  the  lossless  cases,  providing  additional  confidence  of  numerical  model  validity. 

D.  RADIATION  PATTERNS 

Radiation  patterns  were  calculated  for  both  antennas  to  provide  a  final  certification 
of  model  validity.  Elevation  patterns  of  the  crossed  loop  were  obtained  for  2,  5  and  10 
MHz,  under  the  four  ground  conditions,  free  space,  perfect  ground,  good  ground,  and  poor 
ground.  A  horizontal  radiation  pattern  for  5  MHz  and  perfect  ground  conditions  verified 
the  omnidirectionality  of  the  total  gain  in  this  plane.  The  patterns  are  shown  in 
Appendix  C. 

Radiation  patterns  of  the  vertical  monopole  with  wire  mesh  counterpoise  were 
produced  for  10,  15,  and  20  MHz,  under  the  same  four  ground  conditions.  The  elevation 
patterns  are  shown  in  Appendix  D  and  include  the  vertical,  horizontal,  and  total  power 
gain  curves.  An  azimuth  pattern  at  15  MHz  and  perfect  ground  verified  the  omni¬ 
directionality  of  the  total  gain  of  the  monopole.  This  pattern  is  also  shown  in  Appen¬ 
dix  D. 

From  the  preceding,  the  conclusions  are  that  since  long-range  HF  propagation 
generally  occurs  at  higher  frequencies  and  lower  angles,  the  choice  of  using  vertical 
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Fig.  19.-  High  Band  Array  Average  Gain  Versus  Frequency. 


monopoles  for  the  high  band  array  is  good.  Similarly,  the  use  of  the  vertical  crossed  loops 
is  a  good  choice  for  the  low  band  array,  since  shorter  range  HF  propagation  uses  lower 
frequencies  and  high  elevation  angles. 
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V.  HFDF  SSL  ERROR  ANALYSIS 


A.  DF  MATHEMATICAL  MODEL 

If  the  actual  location  of  the  antennas  in  the  "L"  shaped  arrays  is  included  in  the 
equations  for  angles  of  arrival  from  Chapter  n,  the  expressions  can  be  further  simplified. 

Recognizing  that  all  antennas  of  the  array  are  located  at  the  same  height  above  the 
ground  and  placing  the  reference  antenna  which  is  located  at  the  vertex  of  the  array  at  the 

origin  of  the  coordinate  system,  z.  ,  and  become  zero  and  can  be  eliminated 

from  equation  33,  reducing  the  expression  for  elevation  angle  of  arrival  to 


and  the  azimuth  angle  to 


(52) 


These  expressions  can  be  further  simplified  by  aligning  the  first  arm  of  the 
interferometer  with  the  y-axis  and  the  second  arm  with  the  x-axis. 


The  antenna  coordinates  x^  and  become  equal  to  zero,  which  when  applied  to 


variables  defined  previously  result  in 


y^  and  c 


x,y,. 


The  elevation  angle  of  arrival  of  the  wave  can  thus  be  expressed  as 


Since  the  interferometer  compares  the  phase  measured  at  the  antennas  two  at  a  time, 
one  on  each  arm  of  the  interferometer  array,  the  antennas  on  each  arm  are  selected  to  be 
equidistant  from  the  reference  antenna.  Incorporating  these  facts  into  the  analysis  and 
using  the  specific  geometry  of  the  two  antenna  arrays,  the  x-coordinate  of  the  antenna  on 
the  first  arm  of  the  interferometer  is  equal  to  the  y-coordinate  of  the  corresponding 
antenna  on  the  second  arm,  so  jTj  =  y^. 
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Replacing  these  values  in  equations  S3  and  54  and  simplifying,  the  elevation  angle 


of  arrival  becomes 


sin'' 


(  \ 

2 

( 

- 

LiJ 

CO. 

(55) 


or 


V 


cos"' 


'  /o’  ^  «>/ 


(56) 


and  the  azimuth  angle, 


‘t) 


tan'' 


% 


(57) 


B.  BEARING  AMBIGUITY  RESOLUTION 

In  both  of  the  antenna  arrays  the  distance  between  the  reference  antenna  and  the 
first  antenna  of  each  arm  is  less  than  one-half  wavelength  at  the  lowest  frequency  (longest 
wavelength),  so  the  phase  difference  between  these  antennas  will  be  less  than  180°.  As 
explained  in  Chapter  11,  the  phase  difference  is  measured  as  modulo  180°,  therefore  no 
ambiguity  in  the  measurement  will  occur.  A  different  situation  occurs  with  the  rest  of  the 
antenna  pairs  of  the  array  where  the  phase  difference  between  each  element  and  the 
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reference  antenna  can  be  more  than  180°  and  ambiguity  in  the  measurement  of  the  angle 
of  arrival  occurs.  A  means  of  determining  the  number  of  180°  differences  is  needed  for 
finding  the  exact  baseline  phase. 

Each  arm  of  the  interferometer  has  three  baselines  formed  by  the  reference  antenna 
and  each  of  the  antennas  of  the  arm.  The  "short-baseline”  is  formed  by  the  first  antenna 
of  the  arm;  the  "intermediate-baseline"  is  formed  by  the  second  antenna;  and  the  third 
antenna  forms  the  "long-baseline".  To  find  the  exact  phase  difference  of  the  intermediate 
baseline,  the  short-baseline  phase  must  be  measured  first.  The  number  of  180°  phase  rota¬ 
tions  are  computed  using  the  formula 

N  = 

where 

=  short-baseline  phase, 
a.  =  intermediate-baseline  phase  (unresolved), 

=  short-baseline  length, 
d.  =  intermediate-baseline  length,  and 
[  ]  =  nearest  integer  function. 

The  resolved  intermediate-baseline  phase  can  then  be  calculated  by 

<D.  =  yv  (  180°)  +  a  .  (59) 

I  '  '  I 


a 


" 

-  a. 

d 

1 

1  *  ) 

- 

(58) 
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Using  the  same  pioceduie  the  lesolved  long-baseline  phase  can  be  calculated  by 
replacing  intermediate-baseline  values  with  the  corresponding  long-baseline  values. 

Considering  the  validity  of  the  phase  measurements  completes  the  mathematical 
model.  When  the  received  signal  is  not  a  single  plane  wave,  phase  measurements  may 
correspond  to  phases  of  different  components  of  the  incoming  wave  (different  propagation 
modes)  or  to  the  vector  sum  of  several  components  arriving  at  the  antenna.  If  the  resulting 
phase  measurements  differ  greatly  from  one  another,  as  when  two  different  propagation 
modes  are  being  received,  they  do  not  provide  required  information  to  unambiguously 
define  the  angle  of  arrival. 

To  solve  this  ambiguity,  the  system  uses  a  Quasi-Single  Phase  Wave  (QSPW)  test 
to  categorize  measured  phases  into  "phase-linear"  and  "non-phase-linear"  values.  The 
QSPW  test  interpolates  the  resolved  long-baselme  phase  to  the  intermediate-baseline 
phase  and  checks  for  agreement  within  an  acceptable  phase  window.  This  acceptable 
window  is  usually  set  to  30°.  The  computation  performed  on  each  arm  of  the  interfer¬ 


ometer  can  be  expressed  as 


where 


(60) 


O;  =  resolved  Intermediate  Phase, 


=  resolved  long  phase, 


d.  =  intennediate  baseline  length, 


dj^  -  long  baseline  length,  and 


=  acceptable  phase  window. 

If  this  condition  is  satisfied  on  both  legs  of  the  array  simultaneously,  the  azimuth 
and  elevation  calculated  with  these  values  is  said  to  be  a  "phase  linear  frame".  The  wave 
can  then  be  considered  a  Quasi-Single  plane  wave,  which  is  a  good  approximation  of  a 
plane  wave. 

C.  SELECTION  OF  ANGLES  OF  ARRIVAL  FOR  MODEL  COMPARISON 


A  set  of  angles  of  arrival  of  incoming  waves  must  be  defined  to  effectively  test  both 
antenna  arrays.  For  the  azimuth  plane,  a  set  of  angles  that  exploits  die  180**  symmetry  of 
both  arrays  is  selected.  The  resultant  set  of  angles  is  shown  in  Fig.  20. 
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Next,  the  elevation  plane  which  corresponds  to  0  -  90*  in  spherical  coordinates 
is  considered.  Elevation  angles  of  arrival  for  testing  the  arrays  are  suggested  by  the 
tactical  scenario  described  in  Chapter  I  and  the  typical  HFDF  SSL  paths  outlined  in 
Chapter  n.  From  these  conditions  a  set  of  angles  is  chosen  which  thoroughly  tests  the 
array  in  the  elevation  plane,  as  is  shown  in  Fig.  21. 


odo* 


Fig.  21.*  Test  Elevation  Angles 

D.  SSL  ARRAY  PHASE  ERROR  ANALYSIS 

To  find  the  amount  of  error  induced  by  the  pairs  of  antennas  in  the  arrays  in  the 
phase  measurements,  the  selected  azimuth  and  elevation  angles  were  combined  to  define 
directions  of  arrival  of  plane  waves.  These  "test"  waves  are  "incident  waves"  for  NEC 
calculations. 

The  phase  of  currents  induced  in  the  different  antennas  of  the  arrays,  when  excited 
by  each  of  the  incident  waves,  were  calculated  via  NEC  and  compared  with  phase  values 
from  the  mathematical  model.  The  mathematical  model  does  not  include  the  effects  of 
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antenna  coupling,  wave  scattering,  or  ground  conditions,  while  NEC  does  include  them. 
The  difference  between  the  phases  produced  by  the  two  models  is  referred  to  as  the 
"phase  error". 

To  select  a  reasonable  number  of  multiple  combinations  of  azimuth,  elevation, 
frequency,  antenna  pair,  and  ground  conditions,  die  following  cases  were  selected: 

LOW  BAND  ARRAY 

•  Frequency:  5  MHz 

•  Ground  condition:  Good  ground  and  poor  ground 

•  Antenna  pair:  Longest  baseline 

•  Azimuth:  -135°,  -090°,  -045°,  000°  and  045° 

•  Elevation  (90-0):  025°,  035°,  045°,  055°,  and  065° 

HIGH  BAND  ARRAY 

•  Frequency:  15  MHz 

•  Ground  condition:  Perfect  ground 

•  Antenna  pair:  Longest  baseline 

•  Azimuth:  -135°,  -090°,  -045°,  000°  and  045° 

•  Elevation  (90-0):  025°,  035°,  045°,  055°,  and  065° 

The  different  runs  are  summarized  in  die  following  subsec-tions.  Figures  showing 
results  of  phase  error  calculations  are  contained  in  Appendix  E. 
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1.  Low  Band  Array  Phase  Errors 

Fig.  E-1  presents  phase  errors  induced  by  the  longest  baselines  of  arm-1  (y 
axis),  and  arm-2  (x  axis),  of  the  low  band  array,  at  5  MHz  over  good  ground,  for  an 
azimuth  000°  and  varying  elevation  from  025°  to  065°. 

The  elevation  angle  of  the  incident  wave  was  set  to  045°  and  the  azimuth 
angle  varied  between  -135°  and  045°.  Fig.  E-2  shows  the  phase  errors  induced  by  the 
longest  baselines  of  both  arms  of  the  low  band  array. 

Ground  conditions  for  the  above  cases  were  changed  to  poor  ground  and  the 
analysis  repeated.  Fig.  E-3  and  Fig.  E-4  show  the  resultant  phase  errors  for  the  same 
array. 

From  these  figures,  it  is  seen  that  for  the  low  band  array  phase  errors 
calculated  for  atm-1  and  for  arm-2  of  the  interferometer  with  good  ground  conditions 
were  comparable  with  poor  ground  conditions.  This  indicates  that  phase  errors  induced 
by  the  low  band  antenna  array  did  not  change  significantly  with  ground  conditions. 

There  is  a  similar  phase  error  in  both  arms  of  the  interferometer  for  an  incident 
wave  as  an  azimuth  angle  of  045°  as  well  as  from  -135°.  The  symmetry  of  the  incident 
wave,  with  respect  to  both  arms  of  the  interferometer,  is  maintained  for  the  low  band 
elements  the  same  as  it  is  for  high  band  at  these  azimuth  angles. 

The  smallest  phase  errors  occur  at  the  highest  incident  wave  elevation  angles. 
This  conforms  to  the  electromagnetic  properties  of  the  crossed  loop  antennas  which  have 
good  gain  at  high  elevation  angles,  and  poor  gain  at  lower  angles. 
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2.  High  Band  Array  Phase  Errors 

Fig.  E-5  shows  phase  errors  produced  by  the  longest  baselines  of  atm-1  (y- 
axis),  and  arm-2  (x-axis)  of  the  high  band  antenna  at  15  MHz  over  good  ground  and  for 
an  incident  wave  azimuth  of  000®,  and  elevation  angle  varying  between  025®  and  065®. 
Fig.  E-6  shows  phase  errors  from  the  same  baselines  of  the  high  band  array  for  identical 
ground  and  frequency  conditions,  an  elevation  angle  of  045®  and  azimuth  angles  from  - 
135®  to  045®.  The  effects  of  changing  the  ground  conditions  to  poor  ground  are  shown 
in  Fig.  E-7  and  Fig.  E-8.  These  values  for  arm-1  varied  significantly  from  those  calculated 
for  poor  ground  case.  A  similar  effect  is  observed  for  arm-2.  Ground  conditions  do  have 
a  significant  effect  on  phase  errors  for  the  high  band  antenna  array. 

When  the  incident  wave  arrives  from  an  azimuth  of  045®  and  -135®  a  similar 
value  of  phase  error  occurs  in  both  arms  of  the  high  band  interferometer.  This  confirms 
symmetiy  of  the  incident  wave  with  respect  to  arms  of  the  interferometer. 

For  good  ground  conditions  the  phase  errors  of  both  high  band  array  arms  are 
larger  at  high  incident  wave  elevation  angles.  For  poor  ground  conditions  the  largest 
phase  errors  were  at  lower  incident  wave  elevation  angles.  These  results  suggest  that 
stronger  antenna  element  coupling  and  inter-element  scattering  occurs  for  good  ground 
than  for  poor  ground. 

E.  EFFECT  OF  SSL  PHASE  ERROR  IN  AOA  MEASUREMENT 

The  preceding  antenna  phase  errors  do  not  directly  represent  the  final  angle  of 
arrival  (AOA)  errors.  The  actual  azimuth  and  elevation  errors  also  depend  on  specific 
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algorithm  sensitivity  to  these  phase  errors.  As  explained  in  Chapter  n,  a  longer  baseline 
decreases  the  effects  of  phase  measurement  errors  in  determining  the  angle-of-arrival 
(AOA)  For  this  reason  only  the  longest  baselines  were  used  to  determine  the  AOA 
measurements,  with  the  shorter  baselines  used  to  eliminate  ambiguities.  Thus,  the  system 
implements  a  wide  aperUire  array  whfle  maintaining  required  ambiguity  resolution 
capabilities. 

The  magnitude  of  the  AOA  errors  were  determined  by  comparing  NEC  phase  values 
with  those  of  the  math  model.  The  figures  showing  the  results  of  the  AOA  error 
calculations,  for  both  the  low  band  and  high  band  interferometer  array,  are  contained  in 
Appendix  F. 

Fig.  F-1  shows  the  azimuth  and  elevation  error  induced  in  the  HFDF  SSL  system 
by  the  low  band  antenna  array.  These  errors  were  calculated  for  a  frequency  of  5  MHz, 
good  ground  045®  azimuth  incidence  varying  elevation  angles  between  025®  and  065®. 
Fig.  F-2  shows  the  HFDF  SSL  error  in  azimuth  and  elevation  for  the  same  conditions  but 
fixing  the  incident  elevation  angle  at  045®  and  varying  azimuth  between  -135®  and  045®. 
Similarly,  Fig.  F-3  and  Fig.  F-4  show  the  corresponding  azimuth  and  elevation  error  for 
the  low  band  array  over  poor  ground  and  the  same  conditions  of  incident  waves  for  the 
previous  two  cases. 

From  an  analysis  of  these  figures  it  is  observed  that  the  HFDF  SSL  azimuth  and 
elevation  errors  induced  by  the  low  band  antenna  array  are  approximately  one  order  of 
magnitude  less  than  the  individual  phase  errors  from  which  they  are  calculated.  The 
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largest  azimuth  and  elevation  error  of  the  low  band  array  was  less  than  0.2**  for  all  cases 
analyzed. 

Fig.  F-5  summarizes  errors  for  the  high  band  antenna  array  at  15  MHz,  good 
ground,  and  an  azimuth  incidence  angle  of  045**  and  elevation  angles  between  025**  and 
065°.  Corresponding  azimuth  and  elevation  errors  for  a  fixed  elevation  of  angle  of  the  of 
045°  and  varying  azimuth  angles  from  -135°  to  045°  are  shown  in  Fig.  F-6.  Equivalent 
error  calculations  were  made  for  the  high  band  array,  for  poor  ground  conditions  and 
appear  in  Fig.  F-7  and  Fig.  F-8. 

There  is  a  tendency  of  the  system  to  produce  larger  azimuth  errors  at  high  elevation 
angles  of  arrival  of  the  incident  wave  for  the  high  band  antenna  array,  under  good  and 
poor  ground  conditions.  Azimuth  errors  of  die  high  band  array  were  approximately  one 
order  of  magnitude  less  than  phase  errors  for  all  ground  conditions  and  angles  of  arrival 
of  the  incident  wave.  The  largest  azimuth  and  elevation  error  induced  by  the  high  band 
array  was  found  to  be  less  than  0.5°  for  all  cases  analyzed. 

The  lowest  phase  errors  occur  when  the  incident  wave  arrives  at  045°  or  at  -135° 
for  both  the  low  and  high  band  arrays.  This  result  shows  that  the  "L"  shape  geometry  of 
the  antenna  arrays  has  two  preferred  angles  of  arrival  at  045°  and  -135°  for  the  specific 
orientation  of  the  antenna  array  shown  in  Fig.  1. 

In  general,  azimuth  and  elevation  errors  were  slightly  smaller  for  the  low  band 
antenna  than  for  the  high  band. 
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F.  HFDF  SSL  ERROR  AND  SSL  FOOTPRINT 


The  accuracy  of  an  HFDF  SSL  system  is  measured  by  comparing  the  predicted 
transmitter  location  to  the  actual  position.  Emitter  position  errors  can  be  observed  after 
the  application  of  the  vertical  triangulation  technique  to  the  angle  of  arrival  of  the 
incoming  wave. 

Many  factors  can  generate  errors  in  HFDF  SSL  systems.  This  study  was  concerned 
only  with  errors  introduced  by  the  antenna  array;  therefore,  the  lateral  and  range  errors 
introduced  in  the  position  of  the  emitter  by  bearing  errors  found  in  the  previous  sections 
are  considered. 

When  the  virtual  height  of  the  ionosphere  is  assumed  to  be  known  and  in  addition, 
a  no-tilt  case  is  assumed,  the  lateral  error  is  only  a  function  of  the  azimuth  error  intro¬ 
duced  in  the  HFDF  SSL  by  die  antenna  array.  Likewise,  range  error  is  dependent  upon 
the  elevation  angle  error.  Furthermore,  the  lateral  and  range  errors  form  what  is  called  the 
HFDF  SSL  "footprint". 

This  footprint  is  usually  elliptical  in  shape  with  the  size  of  its  major  axis  detennined 
by  range  error,  while  the  minor  axis  is  proportional  to  the  lateral  error. 

The  HFDF  SSL  errors  for  the  system  under  consideration  can  then  be  estimated  by 
applying  the  flat-earth  flat-ionosphere  model  described  in  Chapter  n,  to  the  maximum 
azimuth  and  elevation  errors  calculated  in  the  previous  sectior  for  each  antenna  array. 

A  value  of  0.2°  was  used  for  the  low  band  array,  and  0.5°  for  the  high  band  array. 
Lateral  and  range  errors  of  the  HFDF  SSL  system  were  calculated  for  the  largest  and 


69 


smallest  emitter  distances  of  the  tactical  scenario  and  the  maximum  dimensions  of  the 


minor  and  major  axes  of  the  HFDF  SSL  foo^rint  were  then  determined. 

For  the  low  band  array  the  maximum  dimensions  of  the  major  (range  error)  and 


minor  (lateral  error)  axes  of  the  HFDF  SSL  footprint  ate: 

EMITTER-SSL 

DISTANCE 

PROPAGATION 

MODE 

MAJOR  AXIS 
(RANGE  ERROR) 

MINOR  AXIS 
(LATERAL  ERROR) 

300  km 

E 

4.3  km 

2.1  km 

300  km 

F2 

5.3  km 

2.1  km 

500  km 

E 

8.9  ]un 

3.5  km 

500  km 

7.1  km 

3.5  km  . 

For  the  high  band  array  the  maximum  dimensions  of  the  major  (range  error)  and 
minor  (lateral  error)  axes  of  the  footprint  are: 


EMITTER-SSL 

DISTANCE 

PROPAGATION 

MODE 

MAJOR  AXIS 
(RANGE  ERROR) 

MINOR  AXIS 
(LATERAL  ERROR) 

300  km 

E 

10.7  km 

5.2  km 

300  km 

Fz 

13.1  km 

5.2  km 

500  km 

E 

22.4  km 

8.7  Jun 

500  km 

Fz 

17.7  km 

8.7  km  . 

The  maximum  HFDF  SSL  lateral  and  range  errors  occiur  for  the  high  band  array 
where  the  targeted  emitter  is  located  at  the  larger  distance  specified  by  the  tactical 
scenario,  and  the  incoming  wave  has  an  E-propagation  mode.  The  smallest  lateral  and 
range  errors  occur  when  the  system  uses  the  low  band  array,  the  targeted  emitter  is 
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located  at  the  shortest  distance  specified  by  the  tactical  scenario  and  the  incoming  wave 
has  an  E-propagation  mode. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  characteristics  of  two  HF  antenna  arrays  which  are  part  of  a  phase-interferome¬ 
ter  HFDF  SSL  system  that  operates  in  the  2  to  30  MHz  band  were  analyzed  using  the 
Numerical  Electromagnetics  Code  (NEC)  and  comparing  the  results  with  the  ideal  case. 

The  larger  the  aperture  of  the  antenna  array  compared  with  the  wavelength  of  the 
incoming  signal,  the  better  its  resolution  when  determining  die  angle  of  arrival  of  an 
elevated,  incoming  wave.  Therefore,  in  the  HFDF  SSL  system  of  this  study,  the  longest 
baselines  of  the  phase-interferometer  were  selected  to  determine  AOA  while  the  shorter 
baselines  were  used  to  eliminate  ambiguities. 

The  analysis  included  a  comparison  between  the  effects  of  the  diOTerent  ground 
conditions.  Special  consideration  was  given  to  good  ground  (E^  =  30  and  a  =  0.01 )  and 

poor  ground  (e^  -  5  and  a  =  0.001)  cases.  The  high  band  antenna  array  (10  to  20 

MHz),  suffered  the  maximum  coupling  and  strongest  inter-element  scattering  for  good 
ground,  where  it  produced  the  largest  phase  errors.  For  the  low  band  antenna  array  (2  to 
10  MHz),  changing  ground  conditions  do  not  affect  the  phase  errors  substantially. 

There  is  a  consistently  lower  phase  error  induced  by  the  low  band  array  for  2  to  10 
MHz  than  for  the  high  band  array  using  vertical  monopoles  with  a  wire  grid  mesh 
counterpoise  for  the  10  to  20  MHz  band. 
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The  "L"  shape  of  the  geometry  of  the  two  antenna  airays  had  two  preferred  azimuth 
angles  of  arrival  for  incoming  waves.  These  azimuth  angles  were  045**  and  ‘ISS**, 
respectively  (see  Fig.  1).  Waves  arriving  from  these  angles  induced  an  almost  negligible 
azimuth  enor  in  the  DF  system. 

Angle  of  arrival  errors  for  the  high  band  array  were  less  than  0.5°  and  were  under 
0.2°  for  the  low  band  antenna  system. 

The  dimensions  of  the  minor  and  major  axes  of  the  HFDF  SSL  footprint  were 
calculated  for  the  minimum  and  maximum  emitter  distances  considered  in  the  tactical 
scenario.  The  major  axis  was  called  range  error  and  the  minor  axis  was  called  lateral 
error.  The  maximum  HFDF  SSL  lateral  and  range  error  of  this  system  are  8.7  km  and 
22.4  km  respectively.  These  error  values  occur  when  the  system  uses  the  high  band  array, 
the  targeted  emitter  is  located  at  the  largest  distance  specified  by  the  tactical  scenario  (500 
km),  and  the  incoming  wave  has  an  E-propagation  mode. 

The  smallest  lateral  and  range  error  were  found  to  be  2.1  km  and  4.3  km 
respectively.  They  arise  when  the  system  uses  the  low  band  array,  the  targeted  emitter  is 
located  at  the  shortest  distance  specified  by  the  tactical  scenario  (300  km),  and  the 
incoming  wave  has  an  E-propagation  mode. 

The  Numerical  Electromagnetics  Code  (NEC)  proved  to  be  an  excellent  tool  for  the 
analysis  of  this  electromagnetic  problem.  It  was  particularly  useful  because  of  its  ability 
to  simulate  the  antenna  structures,  the  environment  where  they  are  located,  and  the 
incoming  plane  waves  used  to  test  the  HFDF  SSL  system. 
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B.  RECOMMENDATIONS 


Future  analysis  of  the  effects  of  the  HFDF  SSL  azimuth  and  elevation  errors  found 
in  this  study,  when  using  actual  ionospheric  characteristics  of  the  near-equatorial  region, 
is  needed. 

Comparison  between  theoretical  calculations  performed  in  this  study  and 
measurements  of  the  actual  system  can  provide  information  for  assessing  the  errors 
induced  by  other  components  of  the  system  and  better  insight  into  capabilities  and 
limitations. 


74 


APPENDIX  A 


NEC  DATASETS 


1)  fflGH  BAND  ANTENNA  ARRAY 


CM  HIGH  BAND  7  VERTICAL  MONOPOLES  ARRAY 
CM 

CM  6-rOOT  VERTICAL  MONOPOLES 
CM 

CM  0.5  FEET  ABOVE  THE  GROUND 
CM 

CM  FEED  POINT  AT  BOTTOM  OF  EACH  MONOPOLE 
CM 

CM  GROUND  SCREEN  6  FEET  X  6  FEET  WITH  6  INCH  SPACING 
CM 

CM  ADDITIONAL  6  FEET  RADIALS  EVERY  45  DEGREES 
CM 

CM  SCREEN  UP  0.001  METERS  FROM  GROUND 
CM 

CM  OVER  GOOD  GROUND 


GW111,1, .1524,0. ,0., .1524, .9144,0. ,0.01,  BUILD  ONE  SECTOR 
GW111,1, .3048,0. ,0., .3048, .9144,0. ,0.01,  OF  THE  6  FOOT 
GW111,1, .4572,0. ,0., .4572, .9144,0. ,0.01,  SQUARE  WIRE  MESH, 
GW111,1, .6096,0, ,0., .6096, .9144,0. ,0.01,  0.5  FEET  BETWEEN 

GW111,1, .7620,0.,0., .7620, .9144,0. ,0.01,  WIRES. 

GW111,1, . 9144, 0.,0., .9144,  .9144,0., 0.01, 

GW111,1,0., .1524,0., .9144, .1524,0. ,0.01, 

GW111,1,0., .3048,0., .9144, .3048,0., 0.01, 

GWlll, 1,0., .4572,0., .9144, .4572,0. ,0.01, 

GW111,1,0., .6096,  0., .9144,  .6096,0., 0.01, 

GWlll, 1,0., .7620,0., .9144, . 7620, 0 .,  0.01, 

GWlll, 1,0., .9144,0., .9144,  .9144,0., 0.01, 

GW111,1, .9144, .9144,0. ,1.293, 1.293,0. ,0.01,  CREATE  CORNER  RADIAL 


GMO,3,0.,0.,90.,0.,0.,0.,111.111, 
GW111,4,0. ,-1.8288, 0.,0., 1.8288,0.,  0.01, 
GWlll, 4, -1. 8288, 0.,0.,1. 8288, 0.,0.,  0.01, 
GMO,0,0.,0.,0.,0.,0., 0.001, 111. Ill, 
GWlll, 3, 0.,0.,0. 15, 0.,0., 1.9788, .05, 
GWlll, 1, 0,, 0.,0. 001, 0.,0., 0.15, 0.05, 
GM0,0,0.,0.,0.,0.,4.33,0.,111.111, 
GM1,1,0.,0,,0.,0.,29.0,0.,111.111, 
GM1,1,0.,0.,0.,0.,16.67,0.,112.112, 
GM0,0,0.,  0.,45.,0.,0.,0.,111.113, 
GX3,1C0, 

GM0,0,0.,0.,-45.,0.,0.,0.,111.116, 

GP 

6E1 

GN2, 0,0, 0,30. ,0.01 
FRO, 1,0,  0,15., 

WG, 

NX, 

CM  START  NEXT  STRUCTURE 
CM 

CM  ADD  ANTENNA  ZERO 


GENERATE  THE  OTHER  3  SECTORS 
ADD  Y  AXIS  ADDITIONAL  RADIALS 
ADD  X  AXIS  ADDITIONAL  RADIALS 
MOVE  WIRE  MESH  UP  .001  METERS 
CREATE  6  FOOT  MONOPOLB 
CONNECT  MONOPOLE  TO  WIRE  MESH 
MOVE  STRUCTURE  TO  Y-4.33  M 
CREATE  ANTENNA  2  AT  Y-33.33  M 
CREATE  ANTENNA  3  AT  Y-50  M 
ROTATE  STRUCTURE  45  DBG  UP 
REFLECT  STRUCTURE  OVER  Y  AXIS 
ROTATE  STRUCTURE  45  DEG  DOWN 
SUPPFJSSS  GEOMETRY  PRINTING 
END  GEOMETRY  DEFINITION 
GOOD  GROUND  CONDITIONS 
FREQUENCY  15  MHZ 
CALCULATE  GREEN'S  FUNCTION 
START  NEXT  STRUCTURE 
COMMENTS 
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GW120,1, .1524,0.,0. , .1524, . 9144, 0 . , 0 . 01, 
GW120,1, .3048,0.,0., .3048, . 9144, 0 . , 0 . 01, 
GW120,1, .4572,0.,0., .4572, . 9144, 0 . , 0 . 01, 
GW120,1, .6096,0.,0., .6096, . 9144, 0 . , 0 . 01, 
GW120,1, .7620,0.,0., .762f , . 9144, 0 . , 0 . 01, 
GW120,1, .9144,0.,0., .9144, . 9144, 0 . , 0 . 01, 
GW120,1,0., .1524,0., .9144, . 1524, 0 . , 0 . 01, 
GW120,1,0., .3048,0., .9144, .3048, 0 . , 0 . 01, 
GW120,1,0., .4572,0., .9144, . 4572, 0 . , 0 .01, 
GW120,1,0., .6096,0., .9144, . 6096, 0 ., 0.01, 
GW120,1,0., .7620,0. , .9144,  . 7620, 0 . , 0 .01, 
GW120,1,0., .9144,0., .9144, .9144,0  ,0.01, 
GW120,1, .9144, . 9144,0., 1.293,1. 293, 0.,0. 
GMO, 3,0.,0.,90.,0.,0.,0.,  120 . 120, 

GW120, 4, 0.,-l. 8208, 0.,0.. 1.8288,0..  0.01, 
GW120, 4,-1. 8288, 0.,0.,1. .288,  0 . , 0 . , 0 . 01, 
GMO, 0,0.,0.,0.,0.,0.,0. 001,120.120, 
GW120,3,0.,0.,0.15,0.,0.,1.9788,  .05, 
GW120,l,0.,0.,0.001,0.,0.,q.l5,0.05, 

GP 

GEl 

PTl, 111, 61,0 

EXl, 5, 5, 0,25, -135, 0,10,  45 
XO, 

PTl, 112, 61, 0 

EXl, 5, 5, 0,25, -135, 0,10, 45 


BUILD  ONE  SECTOR 
OF  THE  6  FOOT 
SQUARE  WIRE  MESH, 
0.5  FEET  BETWEEN 
WIRES. 


01,  CREATE  CORNER  RADIAL 

GENERATE  THE  OTHER  3  SECTORS 
ADD  Y  AXIS  ADDITIONAI.  RADIALS 
ADD  X  AXIS  ADDITIONAL  RADIALS 
MOVE  WIRE  MESH  UP  .001  METERS 
CREATE  6  FOOT  MONOPOLE 
CONNECT  MONOPOLE  TO  WIRE  14ESH 
SUPPRESS  GEOMETRY  PRINTING 
END  GEOI4ETRY  DEFINITION 
ANTENNA  1 


ANTENNA  2 


XQ, 

PTl, 113, 61, 0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 114, 61, 0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 115, 61,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 116, 61,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 120, 61, 0 

EXl, 5, 5, 0,25,-135, 0, 10, 45 
XQ, 

EN 


ANTENNA  3 

ANTENNA  4 

ANTENNA  5 

ANTENNA  6 

ANTENNA  0 

END  DATASET 


2)  LOW  BAND  ANTENNA  ARRAY 

CM  LOW  BAND  7  CROSSED  SQUARE  LOOP  ARRAY 
CM 

CM  5 -FOOT  SQUARE  LOOPS 
CM 

CM  0.5  FEET  ABOVE  THE  GROUND 
CM 

CM  ORTHOGONAL  LOOP  2  INCHES  ABOVE  AXIAL  LOOP 
CM 

CM  FEED  POINT  AT  TOP  OF  EACH  LOOP 
CM 

CM  OVER  GOOD  GROUND 
CE 

GWll, 3,0. ,-.762, 1.524,0., .762,1.524, .05,  CREATE  Y  AXIS  LOOP  TOP 
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6W11,1,0., .762,1.524,0., .762,0., .05, 
GW11,1,0., . 762, 0.,0.,-. 762,0., .05, 

6W11, 1, 0., -. 762, 0.,0.,“. 762, 1.524, .05, 
GWll, 3, -.762,0., 1.574, .762,0., 1.574, .05, 
GW11,1, .762,0., 1.574, . 762, 0 . , 0 . 05, .05, 
GW11,1, .762,0., 0.05, -.762,0., 0.05, .05, 
GWll, 1,-. 762,0., 0.05, -.762,0., 1.574, .05, 
GM0,0,0.,0.,0.,0.,0., .15,11.11, 

GM0,0, 0 . , 0. , 0. , 0. , 13 . , 0. , 11.11, 

GM1,1,0.,0.,0.,0.,87.,0.,11.11, 

GM1,1,0.,0.,0.,0.,50.,0.,12.12, 

GM1,1,0.,0.,0.,13.,-150.,0.,13.13, 

GM1,1,0.,0.,0.,87.,0.,0.,14.14, 

GM1,1,0.,0.,0.,50.,0.,0.,15.15, 

GM4,1,0.,0.,0.,-150.,0.,0.,16.16, 

GEl 

GN2,0, 0,0, 30,0.01 
FRO, 1,0, 0,5., 

RT1,11,2,0 

EXl, 5, 5, 0,25, -135, 0,10,  45 
XQ, 

PT1,11,8,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PT1,12,2,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PT1,12,8,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 13, 2,0 

EXl| 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 13, 8,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 14, 2,0 

EXl,  5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 14, 8,0 

EXl, 5, 5, 0,25, -135, 0,10,  45 
XQ, 

PTl  15  2  0 

EXl" 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 15, 8,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 16, 2,0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 16, 8,  0 

EXl, 5, 5, 0,25, -135, 0,10, 45 
XQ, 

PTl, 20, 2,0 

EXl, 5, 5, 0,25, -135,  0,10, 45 
XQ, 

PTl, 20,  8,0 

EXl, 5, 5, 0,25, -135, 0,10,  45 
XQ, 

EW 


CREATE  y  LOOP  SIDE 

CREATE  Y  LOOP  BOTTOl 

CREATE  Y  LOOP  SIDE 

CREATE  X  LOOP  TOP 

CREATE  X  LOOP  SIDE 

CREATE  X  LOOP  BOTTOM 

CREATE  X  LOOP  SIDE 

MOVE  STRUCTURE  UP  .15  METERS 

MOVE  STRUCTURE  TO  Y-13  M 

CREATE  ANTENNA  2  AT  Y-lOO  M 

CREATE  ANTENNA  3  AT  Y-150  M 

CREATE  ANTENNA  4  AT  X-13  M 

CREATE  ANTENNA  5  AT  X-lOO  M 

CREATE  ANTENNA  6  AT  X-150  M 

CREATE  ANTENNA  ZERO 

END  GEOMETRY  DEFINITION 

GOOD  GROUND  CONDITIONS 

FREQUENCY  5  MHZ 

ANTENNA  1  (X  AXIS  ELEMENT) 


ANTENNA  1  (Y  AXIS  ELEMENT) 


ANTENNA  2  (X  AXIS  ELEMENT) 


ANTENNA  2  (Y  AXIS  ELEMENT) 


ANTENNA  3  (X  AXIS  ELEMENT) 


ANTENNA  3  (Y  AXIS  ELEMENT) 


ANTENNA  4  (X  AXIS  ELEMENT) 


ANTENNA  4  (Y  AXIS  ELEMENT) 


ANTENNA  5  (X  AXIS  ELEMENT) 


ANTENNA  5  (Y  AXIS  ELEMENT) 


ANTENNA  6  (X  AXIS  ELEMENT) 


ANTENNA  6  (Y  AXIS  ELEMENT) 


ANTENNA  0  (X  AXIS  ELEMENT) 


ANTENNA  0  (Y  AXIS  ELEMENT) 


END  NEC  DATASET 
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APPENDIX  B 


ANTENNA  INPUT  IMPEDANCE 
TABLE  B-1 


5FT  SO  CROSSED  LOOP  TOP  FEED 


FREQ. 

IMPEDANCE  LOOP-1 

IMPEDANCE  LOOP-2 

MHZ 

(Y  AXIS) 

(X  AXIS) 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

GOOD 

GROUND  (EPS' 

-30,  SIG-0 

.01) 

2 

0.026 

46 

0.02 

45 

3 

0.20 

70 

0.11 

67 

4 

0.55 

93 

0.35 

91 

5 

0.97 

118 

0.66 

115 

6 

1.6 

144 

1.08 

141 

7 

2.2 

171 

1.59 

167 

8 

3.1 

201 

2.27 

195 

9 

4.4 

232 

3.20 

226 

10 

6.0 

266 

4.44 

258 

POOR 

GROUND  (EPS 

-5,  SIG-0. 

001) 

2 

0.59 

46 

0.30 

45 

3 

0.93 

69 

0.51 

67 

4 

1.22 

92 

0.70 

90 

5 

1.51 

116 

0.90 

114 

6 

1.82 

141 

1.12 

139 

7 

2.14 

168 

1.37 

165 

8 

2.49 

195 

1.73 

193 

9 

3.22 

226 

2.16 

222 

10 

3.94 

259 

2.71 

255 

PERFECT 

GROUND 

2 

0.001 

43 

0.001 

42 

3 

0.003 

64 

0.003 

64 

4 

0.011 

86 

0.011 

86 

5 

0.026 

109 

0.026 

109 

6 

0.055 

133 

0.055 

133 

7 

0.10 

158 

0.10 

158 

8 

0.18 

184 

0.18 

184 

9 

0.30 

212 

0.30 

212 

10 

0.47 

242 

0.47 

242 

FREE  SPACE 

2 

0.000 

43 

0.000 

43 

3 

0.002 

65 

0.002 

65 

4 

0.005 

87 

0.005 

87 

5 

0.013 

110 

0.013 

110 

6 

0.028 

135 

0.028 

135 

7 

0.055 

160 

0.055 

159 

8 

0.10 

186 

0.10 

186 

9 

0.17 

215 

0.17 

215 

10 

0.29 

246 

0.29 

246 
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TABLE  B-2 


6  FT  VERTICAL  MONOPQLE.  6FT  SO  WIRE  MESH 


FREQ.  IMPEDANCE 


MHZ 

REAL 

IMAGINARY 

GOOD 

GROUND  <EPS-30,  SIG-0.01 

10 

12 

-541 

11 

13 

-479 

12 

14 

-426 

13 

16 

-380 

14 

19 

-339 

15 

22 

-302 

16 

26 

-268 

17 

32 

-237 

18 

40 

-209 

19 

49 

-184 

20 

59 

-163 

POOR 

GROUND  (EPS-5,  SIG-0.001 

10 

15 

-574 

11 

14 

-511 

12 

14 

-458 

13 

14 

-412 

14 

14 

-371 

15 

15 

-335 

16 

15 

-303 

17 

16 

-274 

18 

17 

-247 

19 

18 

-222 

20 

18  -199 

PERFECT  GROUND 

10 

3.2 

-553 

11 

3.2 

-493 

12 

4.6 

-443 

13 

5.4 

-399 

14 

6.3 

-361 

15 

7.2 

-327 

16 

8.3 

-297 

17 

9.4 

-270 

18 

10.5 

-245 

19 

11.8 

-222 

20 

13.1  -201 

FREE  SPACE 

10 

1.6 

-656 

11 

2.0 

-586 

12 

2.4 

-526 

13 

2.8 

-475 

14 

3.2 

-431 

15 

3.7 

-391 

16 

4.3 

-356 

17 

4.9 

-324 

18 

5.5 

-295 

19 

6.2 

-268 

20 

6.9 

-243 

79 


ASffSMDIX  C 


CROSSED  LOOP  REDIATZON  PRTTERMS 


HFDF  SSL  5FT  CROSSED  LOOP  OVER  PERFECT  GROUND 
ELEVATION  PATTERN  (2MHZ) 


90 


PATTCRN  CAIN  IN  Dll 
■  ■.....■.HOHIZONTAL 

— . .  VtHTICAL 

'——TOTAL 
ELEVATION  ANCU 


HFDF  SSL  SFT  CROSSED  LOOP  OVER  GOOD  GROUND 
ELEVATION  PATTERN  (2MHZ,  EPS=30.  SIG=O.Ol) 


PATTCKN  CAIN  IN  Oil 

•  ■■■■■■■•kohizontal 

—■‘I  ■■■VtHTICAL 
———TOTAL 
ELEVATION  ANCLE 


Figure  C-2 
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Figure  C-3 


HFDF  SSL  5FT  CROSSED  LOOP  OVER  PERFECT  GROUND 
ELEVATION  PATTERN  (5MH2) 
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HFDr  SSL  SFT  CROSSED  LOOP  OVER  GOOD  GROUND 
ELEVATION  PATTERN  (5MHZ,  EPS=30.  SIG=0.01) 


MTTEKN  CAIN  IN  OBI 

■..■.■■■.horizontal 

-VERTICAL 

TOTAL 

elevation  ancle 


Figure  C-5 
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HFDF  SSL  5FT  CROSSED  LOOP  OVER  POOR  GROUND 
ELEVATION  PATTERN  (5MH2,  EPS=5.  SIG=0.001) 


Figure  C-6 
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HFDF  SSL  5FT  CROSSED  LOOP  OVER  PERFECT  GROUND 
ELEVATION  PAHERN  (10MHZ) 


Figure  C-7 
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HFDF  SSL  SFT  CROSSED  LOOP  OVER  GOOD  GROUND 
ELEVATION  PATTERN  (10MH2,  EPS=30.  S1G=0.01) 


Figure  C-8 
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HFDF  SSL  SFT  CROSSED  LOOP  OVER  POOR  GROUND 
ELEVATION  PATTERN  (lOMHZ,  EPS=5,  SIG=0.001) 
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Figure  C-9 
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HFDF  SSL  SFT  CROSSED  LOOP  OVER  PERFECT  GROUND 
HORIZONTAL  PATTERN  (5MH2) 
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•VWTICAL 

•TOTAL 


ANCLES  IN  SCGRCES  TRUE 


Figure  C-10 
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APPENDIX  D 


VERTICAL  MONOPOLE  RADIATION  PATTEBNS 


HFPr  SSL  6FT  MONOPOLE.  6FT  SQ  WIRE  MESH.  OVER  PERFECT  GROUND 
ELEVATION  PAHERN  (10MHZ) 


90 


180 


Figure  D-1 
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HFDF  SSL  SFT  MQNOPOLE.  SFT  SQ  WIRE  MESH.  OVER  GOOD  GROUND 
ELEVATION  PATTERN  (10MH2,  EPS=30.  SIG=0.01) 


PATTCRN  CAIN  IN  Oil 
.■■■■■.■•.horizontal 
•  ■  ■  VERTICAL 

—————TOTAL 

ELEVATION  ANCLE 


Figure  D-2 
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HFDF  SSL  SFT  MONOPOLE,  6FT  SQ  WIRE  MESH.  OVER  POOR  GROUND 
ELEVATION  PATTERN  (lOMHZ,  EPS=5.  SIG-0.001) 


Figure  D-3 
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HFDF  SSL  6FT  MONOPOLE.  6FT  SQ  WIRE  MESH.  OVER  PERFrCT  GROUND 
ELEVATION  PAHERN  (15MH2) 
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HFDF  SSL  6FT  MONOPOLE.  SFT  SQ  WIRE  MESH.  OVER  POOR  GROUND 
ELEVATION  PATTERN  (15MHZ.  EPS=5,  £10=0.001) 
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Figure  D-6 
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HFDF  SSL  SFT  MONOPOLE.  6FT  SQ  WIRE  MESH.  OVER  GOOD  GROUND 
ELEVATION  PATTERN  (20MHZ.  EPS=30,  SIG=0.01) 
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Figure  D-8 
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HFDF  SSL  6FT  MONQPQLE.  6FT  SQ  WIRE  MESH.  OVER  POOR  GROUND 
ELEVATION  PATTERN  (20MH2,  EPS=5,  510=0.001) 
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Figure  D-9 
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LOW  BAND  ARRAY  PHASB  ERROR  <5  MHz 


Figure  E-2 
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Figure  E-3 
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Figure  E-4 
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Figure  E-5 
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Figure  E-6 
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Figure  E-*7 
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HIGH  BAND  ARRAY  PHASE  ERROR  fl5  MHz 
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Figure  F-2 
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Figure  F-4 
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